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Abstract 

In this work we have developed the essential tools for the algebraic Bethe ansatz solution 
of integrable vertex models invariant by a unique U(l) charge symmetry. The formulation is 
valid for arbitrary statistical weights and respective number N of edge states. We show that 
the fundamental commutation rules between the monodromy matrix elements are derived 
by solving linear systems of equations. This makes possible the construction of the transfer 
matrix eigenstates by means of a new recurrence relation depending on N — 1 distinct types 
of creation fields. The necessary identities to solve the eigenvalue problem are obtained 
exploring the unitarity property and the Yang-Baxter equation satisfied by the i2-matrix. 
The on-shell and off-shell properties of the algebraic Bethe ansatz are explicitly presented in 
terms of the arbitrary i?-matrix elements. This includes the transfer matrix eigenvalues, the 
Bethe ansatz equations and the structure of the vectors not parallel to the eigenstates. 
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1 Introduction 



The quantum inverse scattering method has been playing a major role in the development of the 
theory of two-dimensional integrable models [H[2j[3]. This method paved the way for the discovery 
of important generalizations of the six-vertex model [I] as well as it has helped to prompt the 
notion of quantum group symmetry p2, [6] . This approach also offers us the basic tools to solve 
exactly quantum integrable systems including the computation of lattice correlation functions 

Em. 

The central object in the quantum inverse scattering method is the monodromy matrix T^(X) 
which depends on the continuous spectral parameter A. This operator is frequently viewed as a 
matrix on the auxiliary space A and for an arbitrary ^-dimensional space one can write, 

N 

TAX) = £ r ^( A Kf > C 1 ) 

a,b=l 

where e a ^ are the standard N x N Weyl matrices. 

The matrix elements 7^&(A) act on the space of states of some quantum physical system and 
they are the generators of a quadratic algebra denominated Yang-Baxter algebra. The set of 
relations defining this algebra are, 

N N 

£ R(\,rf$r f MTM = £ T bJ ^)T a , g (X)R(X^) c g j , (2) 
/,<?=! f,g=l 

In analogy to Lie algebras, functions -R(A, /i)^ can be interpreted as the structure constants 
of the Yang-Baxter algebra (T5J). They are often viewed as the elements of a N 2 x iV 2 .R-matrix 
acting on the tensor product of two auxiliary spaces. This matrix can be defined as, 

N 

-Ri2(Ai,A 2 )= £ -Rq^Ai, A 2 )e ajC ® e M . (3) 

a,b,c,d=l 

The associativity of the Yang-Baxter algebra requires that the i?-matrix satisfies the celebrated 
Yang-Baxter equation [H], 

-Ri2(Ai, A 2 )-Ri3(Ai, A 3 )-R 2 3(^2, A 3 ) = -R 2 3(A 2 , A 3 )i?i 3 (Ai, A 3 )i?i 2 (Ai, A 2 ), (4) 
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where R a b(^a, A&) denotes the i?-matrix acting on the tensor product of the spaces A a <8> Ab- 

In this paper we shall consider the solutions of the Yang-Baxter equation (j3J) that can be 
normalized in order to satisfy the unitarity property, 

R 2 i(\,ij i )Ri2(h,\) = In®In, (5) 

where 1^ is the N x N identity matrix. 

At this point we remark that the unitarity property ([5]) follows from the Yang-Baxter equation 
(jlj) under the extra assumption that the i?-matrix is regular. This hypothesis is equivalent to say 
that i? a fe(A, A) is proportional to the permutator on ® . The unitarity assures us that the 
i?-matrix is invertible and this property is important to show that the Yang-Baxter algebra leads 
us to mutually commuting operators. This family of commuting operators is then obtained by 
taking the trace of the monodromy matrix on the auxiliary space, 

N 

T(A) = ^T a , a (A), (6) 

a=l 

In the theory of integrable models the operator T(A) is regarded as the generating function 
of the corresponding quantum integrals of motion such as the underlying one-dimensional Hamil- 
tonian. The understanding of the physical properties of these quantum systems should therefore 
include at least the knowledge of the eigenvalues and the eigenvectors of T(A). In principle, the di- 
agonalization of T(A) can be accomplished through the formulation of the algebraic Bethe ansatz. 
For comprehensive reviews on this subject see for instance the references [31 [TU] - The basic idea of 
this method is to exploit commutation relations between the monodromy matrix elements coming 
from the Yang-Baxter algebra ([2]). In particular, the eigenstates of T(A) are constructed by apply- 
ing appropriate off-diagonal monodromy matrix elements, usually named creation operators, on a 
previously chosen reference state. The main expected feature of this state is that the action of the 
monodromy matrix on it gives us as result a triangular matrix. We emphasize, however, that the 
existence of such reference state does not immediately guarantee the success of an algebraic Bethe 
ansatz solution. In fact, for general values of N, we are not aware of any recipe to implement the 
algebraic Bethe ansatz even when a possible reference state is the trivial ferromagnetic highest 
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vector. For this class of models, our current knowledge on algebraic Bethe ansatz formulations 
remains restrict to very specific integrable systems. A representative category of such models 
are those directly related to i?-matrices based on the vector representation of some special Lie 
algebras [TH [121 13] an d ^2 graded superalgebras [HI [151 116] - 

In order to bring new insights into the algebraic Bethe ansatz approach we have to consider 
the diagonalization of T(A) without referring to any particular dependence of the i?-matrix on the 
spectral parameters. This point of view is rather evident when the i?-matrix commutes with at 
least one U(l) charge symmetry. This feature assures us that the ferromagnetic pseudovacuum 
will play the role of a suitable reference state for the most fundamental integrable model associated 
to such i?-matrix. The diagonalization of the corresponding T(A) should then be completed solely 
on basis of the commutation relations derived from the Yang-Baxter algebra and the constraints 
imposed by both the Yang-Baxter equation and the unitarity property. Despite of its relevance, 
this strategy of solving integrable models has so far remained largely ignored in the literature. 
This is particularly the case when the dimension of the monodromy matrix is iV > 3 since we 
have to consider the presence of different types of creation fields. The basic problem is to unveil 
the role played by each creation operator on the structure of the eigenvectors. This understanding 
is certainly more difficult when the i?-matrices elements are not specified. It is expected that 
identities among the weights will be crucial to solve this problem. 

In this paper we report on some progress towards the algebraic Bethe ansatz solution of 
integrable models with arbitrary .R-matrix []. We will consider the simplest possible family of 
models whose Hilbert space description requires us to consider many independent quasi-particle 
excitations. This turns out to be the systems that their i?-matrices commute with a single U(l) 
symmetry for arbitrary values of N, 

[R 12 {\,fj),S z ® I N + I N ®S Z ] = (7) 

where S z denotes the azimuthal component of an operator with spin s = (N — l)/2. Note that 
this invariance implies that R(X, fj)^t ^ om y wnen the ice rule a + b = c + d is satisfied. 
1 A small part of our results has been briefly announced in reference |17) . 
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We can use property (J7J) to express the i?-matrix in terms of the sectors labeling the 2N — 1 
eigenvalues of the U(l) operator. In the Weyl basis indices these sectors are easily parameterized 
by the charge q = a + b — 1 and the i?-matrix can be written as, 

2JV-1 m{q,N} 
9=1 a,c=M{l,q+l-N} 

where M{x, y} and m{x, y} denotes the maximum and minimum integer of the pair {x, y}, re- 
spectively. 

These integrable systems can be seen as multistate extensions of the totally asymmetric six- 
vertex model preserving a single conserved quantum number. The classical example is the model 
whose i?-matrix is based on the higher spin representation of the quantum U q [SU(2)] algebra 
[T8] . Here we stress that a systematic classification of U(l) invariant solutions of the non-linear 
Yang-Baxter equation is beyond the reach at present. It is therefore rather probable that a variety 
of U(l) invariant i?-matrices are still waiting to be discovered. This should include solutions with 
non-additive /^-matrices as well as those with infinity number of degrees of freedom typical of 
non-compact models. An example of the latter system could be the vertex model based on the 
discrete infinite dimensional representation of the SL(2, R) algebra. Note that such representation 
has a highest weight which plays the role of possible reference state. The framework developed 
here has the nice feature of being able to accommodate the algebraic solution of all such different 
systems in a rather unified way. 

We have organized this paper as follows. In next section we describe the basic properties of 
the vertex model associated to the fundamental representation of the Yang-Baxter algebra based 
on the R- matrices satisfying the ice rule ATI) . In section 3 we describe the essential tools that are 
necessary to obtain the appropriate commutation rules between the monodromy matrix elements. 
It is argued that these relations are obtained by solving coupled systems of linear equations. To 
fix up the main idea of our method we presented the structure of specific sets of commutation 
rules. In section 4 we discuss a procedure to obtain suitable identities between the elements of 
the i?-matrix from that Yang-Baxter and unitarity relations. These identities are decisive for the 
solution of the transfer matrix eigenvalue problem. In section 5 we diagonalize the transfer matrix 
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by means of an algebraic Bethe ansatz. The respective eigenstates are constructed similarly to 
that of a bosonic Fock space with N— 1 distinct creation fields. We present the explicit expressions 
for both the on-shell and off-shell parts of the action of the transfer matrix on the eigenstates in 
terms of the statistical weights. Our conclusions are presented in section 6. In appendices A and 
B we summarize the technical details entering the solution of the two and three particle problems, 
respectively. 

2 The vertex model representation 

Each solution of the Yang-Baxter equation (j3j) gives rise to representations of the Yang-Baxter 
algebra (j2j). The representations depending on the spectral parameter are named Lax operators 
denoted here by C^i(X). The simplest such representation can be obtained directly from the 
i?-matrix amplitudes by the expression, 

N 

C M (X) = Yl fl&jOjjea.eOeg, (9) 

a,b,c,d=l 

where e^ d are N x N Weyl matrices acting on the tensor product Ylf=i ®f space of an one- 
dimensional lattice of length L. The variables /ij play the role of free continuous parameters. 

The Yang-Baxter algebra has the co-multiplication property that the tensor product of two 
representations is still another possible representation. Hence, the following ordered product of 
/^-operators 

Xa{\) = Cal{\ Hl)Cal-i{\ (J>l-i) ■ ■ ■ £ai(A, Hi), (10) 

is indeed a representation of the quadratic algebra (|2j). 

In the context of classical vertex models of statistical mechanics Cj^{X) encoded the structure 
of the local Boltzmann weights at the i-th site of a square lattice of size L. The possible states of 
such statistical system are then associated to the possible bond configurations a,b,c,d— 1, . . . , N 
of each vertex on the Lx L lattice. The energy of each i-th vertex configuration is then associated 
to the statistical weight R(X, Hi) a ,d b - The row-to-row transfer matrix of these vertex models T(A) 
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can be written in a compact form with the help of the monodromy matrix (fit)]) . It is given by the 
trace of T^(A) over the auxiliary space, 

T(A) = Tr A [C AL (\, Hl)Cal-i(K^l-i) ■ ■ ■ Aai(A, Hi)]. (11) 

In order to diagonalize the operator T(A) within the quantum inverse scattering framework 
one needs to know an exact eigenstate |0) of T(A). This vector works a reference state in the 
construction of the Hilbert space of T(A) by an algebraic Bethe ansatz. In general, such state 
is searched by asking that the action of the lower left monodromy matrix elements on |0) are 
annihilated for arbitrary A. This means that the monodromy matrix acts as an upper triangular 
matrix on |0), namely 

w a (A) |0) , for a = b 

Ta,bW |0> = I 0, for a > b (12) 

\ab) , for a < b, 
where \ab) denotes a general non-null vector. 

The action of the upper and the lower elements T a ,b{X) on |0) has therefore distinct meanings. 
For a < b they play the role of creation fields while for a > b they can be thought as annihilators. 
The presence of the U(l) invariance ([7]) makes it possible to build up such reference state in terms 
of the tensor product of local ferromagnetic vectors, 



io>=n®i s >; 



i«>< 



i=i 







V / 



(13) 



JV 



The state \s) i can be viewed as the highest eigenstate of the spin operators 5*f acting on the 



i-th site of a chain of length L. The total spin of the reference state is therefore 



(N-1)L 



and this 



is the reason why |0) is called the ferromagnetic state. From Eqs. (!8ll9][T0l) it is not difficult to see 
that property (TT2"|) is in fact satisfied and that the expressions for functions w a (X) are, 



«>a(A) = U-R(A,//i)£}. 



(14) 



i=i 



The next step would then be the construction of other eigenstates of T(A) besides |0). To this 
end we shall first consider the action of the total spin operator Y^i=i 011 the nionodromy matrix 
elements T a ^ (A). Considering the structure of the Lax operators (|8|9l) it is possible to derive the 
following relation, 



i=l 



(6-a)T 0)6 (A). (15) 



The above commutation relation is useful to illuminate the physical content of the creation 
fields T a ^(X) for a < b. In fact, by acting Eq. ffl5l) on the reference state |0) one derives the 
property, 

L 

Sir a , b {\) |0> = [s(L - 1) + s - b + a] T a , 6 (A) |0) for a < b. (16) 



i=l 



From Eq. (1161) we conclude that the field T a ^{\) for a < b behaves as raising operator of an 
excitation over |0) whose value of the azimuthal spin component is s a> b = s — b + a. This means 
that all monodromy matrix elements of a given upper diagonal have the same azimuthal spin and 
therefore describe the same type of excitations. As a consequence of that we end up with only 
N — 1 linear independent creations fields which is represented, in the simplest manner, by the first 
row of the monodromy matrix Ti^X) for b = 2, • • • , N. The analogy with the Hilbert space of high 
spin Heisenberg magnets strongly suggests that the eigenstates of the transfer matrix (TTTT) could 
be build up in an algebraic fashion by using the help of such N — 1 independent creation fields. 
The whole construction, however, will depend very much on our ability to recast the Yang-Baxter 
algebra (T5]) in the form of convenient commutation rules between the diagonal and the off-diagonal 
monodromy matrix elements. In next section we will deal with this artistic part of the quantum 
inverse scattering method. 



3 The fundamental commutation rules 

The purpose of this section is to discuss the structure of the commutation relations that are 
relevant in the exact diagonalization of the transfer matrix ffTT]) . In principle, these relations 
between the monodromy matrix elements are derived considering the Yang-Baxter algebra ([2]). 
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Let us use the symbol [a; j3] to represent the a-th row and the /3-th column of the iV 2 x iV 2 matrix 
defined by Eq.([2]). The projection of the Yang-Baxter algebra on the rather general set of entries 
[(a — 1)N + b; (c — 1)N + d] can be written as, 

m{a+b-l,N} m{c+d-l,N} 
e=M{l,a+b-N} e=M{l,c+d-N} 

a,b,c,d = 1 . . . N, (17) 

In what follows we will discuss three different families of commutation rules between the 
diagonal, creation and annihilation operators which can be derived from Eq. (flT|) . 

3.1 The diagonal and creation fields 

The creation operators (A) for b = 2, . . . , N will provide us a basis to construct the eigenvectors 
of the transfer matrix T(A). Their commutation rules with the diagonal monodromy matrix 
elements T aA (X) for a = 1, . . . , N are therefore essential in the transfer matrix eigenvalue problem. 
These relations are obtained from the entries [a; (a + c — l)iV + b — c] upon suitable choices of the 
extra variable c. By substituting a = 1, b = a, c = a + c and d = b — c in Eq. lflTl) we find that the 
respective equations are, 

a m{a+b-l,N} 
e=l e=M{l,a+b~N} 

In general, to obtain commutation rules that are useful for the eigenvalue problem, we still have 
to elaborate on Eq. (fl8|) . The additional manipulation consists in making particular combinations 
of a number of equations derived from Eq. (fl8l) with the help of the variable c. In Table 1 we have 
summarized in detail the linear combination required for each diagonal field. 

From Table 1 we see that only the commutation rules for the operators ^^(A) and T^^(X) 
follows directly from Eq. ffl8|) . Indeed, by setting c = b — 1 for 7^ (A) and c = for T N<N (\) one 
finds that their commutation rules with 7^ (/x) are, 

TiAWM = ^lM. r 1 MT ltl (X) - £ ^^^-(A)^^) (19) 

A) M g=2 K(n, A) b ' A 
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Operator 


Diagonal 


Combination 


Number of 




index 


index 


equations 


T hl (n)T ljb (\) 


a = 1 


c = b- 1 


1 


T a!a (X)T ltb (fj,) 


2 < a < N + 1 - b 


c = 0,...,6-l 


b 


T aja (X)T 1)b (fj,) 


N + 2 - b < a < N 


c = 0, . . ., JV- a 


N + l-a 


T N , N (\)T 1:b (n) 


a = N 


c = 


1 



Table 1: The linear combination made from Eq. (fT5)) to obtain commutation rules between the fields T a ^ a {\) and 
Tl^A). The number of equations entering the linear combination is governed by the index c. 



T N , N (\)Ti, b {ii) 



N R(\,nf' h 



— — - w ^l 1>b (n)l NjN (A)+ ^ . . N>1 

K K A iWn,1 e=b+l K \ X iNNA 



\-b _£z£^ g(yLt)7iV ,JV+b-e(A) 



N-l 



N.l rr IWrr I \ 
— aNJ % jN (\)T N - S +l,b(fJ0 



e=l 



(20) 



The commutation rules for the remaining diagonal operators T aA (X) for 2 < a < N — 1 
demand a considerable amount of extra work. In these cases we first have to implement the linear 
combination according to the indices c exhibited in Table 1. As an example let us consider the 
linear system associated to the diagonal fields with 2<a<N + l — b. We find that the relations 
coming from Eq. (fT8]) for c = 0, . . . , b can be arranged in the following form, 



4 a ' b) (A,/i) 



( r x M%^b-x{\) ^ 

^1,2(^)^,0+6-2 (A) 



y T 1>b {ij)T a>a {\) J 

a+b— 1 



$^(A,/i)aJ 



e+l 



e=l 



f 7"g )a (A)T a _e+l,fe(/i) ^ 
%,a+l W%,-e+l,b-l (aO 



y / ^g,a+6-l(A)Zj_g + i ) i(/i) J 



- ^ a ' 6) (A^)T lje »T a)a+b _ e -(A) for 2<a<iV + l-6, (21) 

e=fe+l 
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where the b x b matrix A^' b \\, jj) is given in terms of the structure constants as 



At b \\^) 



^ R{\ ^Td+b-1,1 R (K VYa+b-2,2 
a+6-1,1 



\a+6-l,l o/\ \a+b— 1.1 n / \ \a+&— 1,1 

■ R\^Na,b 



\ R(\ HYa+b-il R {\ f^l+b-2,2 



(22) 



while the 6x1 vector v[ a ' b \\, jj) is 



/ 



(A,//) 



\ 



a+1,6-1 

e,e 



(23) 



We have now reached a point in which the properties of linear systems of equations are appli- 
cable. In fact, the term Ti^{fJ>)%,,a{X) can be systematically single out from the left-hand side of 
Eq. (j2ip with the help of Cramer's rule. Besides that we note the first non-homogeneous part of 
Eq. fTSTj) contains the reverse product 7^ a (A)7^ 6 (//) when the sum index is e = a. Therefore, by 
computing determinants of b x b matrices we are able to generate a unique linear equation relating 
the operator products 7^ a (A)7i i b(/i) and Ti^(ij)T a ^(\) . This equation is then easily solved in order 
to establish the commutation rule between T aA (X) and Ti^(n) for 2<a<N + l — b. 

Clearly, the same method described above can be used to obtain the commutation rules for 
the other values iV + 2 — b < a < N — lof the index a. The system of linear equations is however 
different and for sake of completeness we also present it here. Following Table 1 one finds that 
the respective Eq. ffT8l) can be written as, 

/ %,a{\)%-e+X,b{^) ^ 
Tl.a+b-N+l{^)T~a,N-l{X) v-^ xea-e+1 %,a+l (A)Zi-e+l,b-l (fl) 



1 T 1<a+b _ N (n)T a>N (\) \ 



e+1 



e=l 



y T 1>b (ij,)T ata (X) J 



I %,NW%i-e+l,a+b~-N(^) J 



N 



^°' 6) (A,//)r lie -(//)T aia+6 _e(A) for N + 2-b<a<N-l, (24) 



e=b+l 
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where the (N + 1 — a) x (JV + 1 — a) matrix A2 •'(A, /i) is given by 



i?(A,/x) a '" 



4 a ' b) (A,/,) 



2+6— JV 

i?(A,/x) 

N,a+b-N 



p/\ \N,a+b-N 
\ n \ A J ^)N,a+b-N 



R(X,fi) 



N-l,a+b-N+l 

a+1,6-1 
AT-l,a+6-W+l 



#(A,/x) 



a+1,6-1 

a.b 



T}/\ \N,a+b-N 



d/x \N,a+b-N 



(25) 



while the (iV+l-o) x 1 vector (A, /x) is 



,(«.&) 



(A,/i) 



r>/\ \N,a+b-N 
\R(\fj)a+b-e,e ) 



(26) 



As before the desirable commutation relations between T a>a {X) and ^^(/x) for iV + 2 — b < a < 
N — 1 are obtained by using Cramer's rule in the linear system (I24H26I) . We emphasize that the 
commutation rules discussed in this subsection have the property that the corresponding 7j g (/x) 
fields that appear on the left-hand side of the products satisfy the condition e > b. Here we recall 
that these type of fields are the ones with clear potential to participate directly on the structure of 
the eigenstates. The condition e > b guarantees that such operators, generated by taking 

T aa (X) through 7^(/x), will only contribute to eigenvectors made out creation fields with spin 
lower or equal to s\^. In this manner we assured that the eigenvectors constructed in terms of the 
N — 1 7i fc(/x) operators can indeed be thought as multiparticle states ordered by their spin values. 

Let us now illustrate how the general method explained above works in practice. To this 
end we shall present explicitly the complete structure of the commutation rules between the first 
two creation fields with all the diagonal operators. The simplest example concerns to the 7^2 (/x) 
operator. Its commutation rules with the diagonal operators play a fundamental role in the transfer 
matrix eigenvalue problem since they dictate in which way the transfer matrix eigenvalues depend 
on the arbitrary i?-matrix elements. These relations for a = 1 and a = N follows directly from 
Eqs. (11911201) while those for 2 < a < N — 1 are obtained by solving the 2x2 linear system of 
equations (l21H23p . It turns out that the final results for the commutation rules between 7^ a (A) 
and T lt2 (ij) are, 
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R(fi, X) 1 \ 

.,2^ 1,2(^1,^) 
#(/i,A) 21 



R(fi, A) 2 ' 1 



0+1 



A) 2 \ 



r 7I i2 (A)7l,i(/x), 



(27) 



e=3 



o-l /l) 5,a_e+1 

^ "'^,1 ^,o(A)^_e+l,2(At) 

e=l 



e,a+l 



for 2<o<iV-l, 



(2? 



T NiN {\)T 1>2 (ij) 



R(X,H) 



N,2 
N,2 



jT lt2 (fi)T N)N (X) + 2J 



-R(A, aOam 



>N+2-e,e 



Y^1~i !S (fJ,)TN j N+2-e(X) 



e=3 



„ /x ^jv a 7i )A r(A)7/v_g + l,2(Af) 



e=l 



(29) 



Function _D 2 a ' e ^(A, //) is defined by using the determinant of a matrix whose elements combines 



the first column of A ( f' 2 \x, fi) with the vector v \ a ' z> (X, fi). Its expression in terms of the i?-matrix 
elements is, 



o+l,l -R(A,/i)a^g,e+2 



o.2 



0+1,1 



e = 0,---,a-l and 2<a<iV-l. 



v O+l.l 



for 



(30) 



The next simplest case is related to the T^/x) operator. In order to get the complete set 
of commutation relations we have now to solve two different linear systems with sizes 3x3 and 
2x2. According to Table 1 the former is defined by Eqs. (12 111231 while the latter is associated to 
Eqs. (l24t l 26|) . By solving such linear systems of equations we find that the commutation relations 
between T aa (X) and T 13 (//) are given by, 



R(n, X)ii 

— ^^1,3(^)^1,1 (A) 



R{^i A)a2 



#(^A) 3 \ 



^r 1)2 (A)r li2 (//) 



R(li, X)V 



'^{X)T hl {n) 



(31) 



12 



a+2 



a— 1 n/\ ..\e,a.— e+1 



e=4 



E 

e=l 



^(A,/i):;r 



-Tg ia (A)r a _ e - +li3 (/x) 



+ ^^(A,/i):r +i 



i?(A,/x)^ 21 i?(A,/i)^ 12 



a, 3 



e=l 



-R(A,/x)^i 

\a+2,l 



-R(A, /i)"+};2 

\a+2,l 



R(\, /J,)l+2,1 R{\nTail,2 



a, 3 



a+1,2 
a+1,2 



R{\^1 



e=l 



-R(A,/i)„+2a i?(A,//)^2 



a+1,2 



for 2 < a < iV - 2, 



'7g,a+l(A)Zi-e+l,2(A t ) 



/ 7g,a+2(A)7^_e + l,l(/x) 



T^-^jv-^A)^^) = 



X 



R{\H)N,2 



i?(A,/i)f' z 



R(X,fi) 



N,2 
N-1,3 



N 



- E 



R{X, /i)^_J^i?(A, 

R(^^)n,2 1,3 R(\ V)N+2-e,e 



Ti !3 (//)7^_i iJV -i(A) 



^l,e(A t )^7V-l,Af+2-e(A) 



e=4 
JV-2 



+ E 



e=l 



R(X, fj)^_\'^R(X, n)% ,2 2 
R{X, /^)ty 2 1 ' 3 

— — - N - 2 -%,N{X)T N _ s ,2(fJ') ~ %,N-l{X)T N _^{lj) 

R{\Nn,2 



R{\jj^Zi^ 
R(\^)n-ia 



+ 



R(X,fi) 

^r2—T N ^ ltN (X)T lt2 (n) 



^(A,/f)jv' 2 
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^v,7v(A)Ti i3 (/i) = 



R^iHYnS^- i \rr i\\ , v- />'(A.//) V 

^^i,3(w^,jv(A) + 2^ 



,Af,3 



-R(A, aO^i 



e=4 



, v | Ti i e(/i)7^ )Ar+3 _g(A) 



N-l 



~ n/\ .^AT,l ^Mv(A)7n-S+1,3(aO 



e=l 



-R(A, aOam 



(34) 



where function D^ a ' e \\, fi) is represented by the ratio of 3 x 3 and 2x2 determinants, namely 

R(\,H)l+ 2A R(\,(j)%Zi t2 R (\^)a-e,3+e 

R(\,V) 



for 



#(A,/i)a+ 2 ;i -R(A,^)a^2 -R(A,/i)a +2 



e,3+e 



^(A,/i)aS;i ^(A,/^)atl'2 



^(A,//)a+|l ^(A,/i)a^2 

0, • • • , a - 1 and 2 < a < N - 2 



\a+2,l 



(35) 



3.2 The creation fields 

The results from previous subsection reveal us that the commutation rules between the opera- 
tors T aA (\) and T lib (\) are able to produce additional creation fields other than the basis vectors 
Ti t b(X) or Ti t b(fi). Though these extra creation operators do not take a direct part on the multi- 
particle state basis they are essential for the solution of the transfer matrix eigenvalue problem. 
It is therefore necessary to disentangle not only the commutation rules among the basis vectors 
themselves but also those involving the operators Tl,&(A) with the remaining creation fields T a ^(/i) 
for b > a — 2, • • • , N — 1. In the course of our analysis we find convenient to describe such com- 
mutation relations by using the fields Ti^-^A) and T^-i^+^/i) where the underlying indices 
belong to the following intervals, 

2 < ai < N, 0<di<JV-ai, 2<b = b 1 -d 1 <N (36) 

The above arrangement of indices bring us at least two technical advantages. First, we assure 
that all the mentioned commutations rules between the creation fields will be considered without 
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unnecessary repetition. Next, Eg. ([16]) implies that the effective azimuthal spin component asso- 
ciated to the term Tx^XjT^-i^+d^) should be indexed by the composed variable b + d\. The 
parameterization b = b\ — d\ has therefore the merit of allowing us to keep track of the azimuthal 
spin of such product of creation fields in terms of a unique index bx . As a consequence of that the 
number of the distinct commutation relations that are needed in a given multiparticle state sector 
can now be controlled with the help of the remaining indices a\ and d\. A systematic analysis of 
the Yang-Baxter relation ( IT7T) reveals us that these commutation rules are derived from the entries 
[di — 1; (ai + c — l)iV + bi — c] for selected values of the variable c. The corresponding relations 
are obtained from Eq. (fT71) by choosing a = l,b = ai — 1, c = ai + c e d = bi — c, namely 



oi-l m{ai+6i-l,JV} 

22 R{\^Ta^-ll%Ai+c{^)T'a 1 -e,b 1 -c{p) = ^ 7l j g(^)7^ 1 _l j0l+ 6 1 _e(A)i2(A, /x)^+bi-g,g- 

e=l e=M{l,a 1 +b 1 -N} 

(37) 

The majority of the suitable commutation rules between the fields Ti^-d^A) and 
T ai -x jai+ di (/-*) require a considerable amount of manipulations among the relations (1371) . As before 
we still have to implement certain linear combinations by using the freedom of the index c. This 
procedure is highly dependent on the variables ai, b\ and it has been detailed in Table 2. 

The simplest type of commutation rules are those involving the creation fields that has direct 
participation on the eigenvector basis. These relations are sorted out by choosing the index a% = 2. 
From Table 2 we see that the commutation rule between 7i ) j 1 _d 1 (A) and Ti^+^if^) for b\ > N 
follows directly from the entry [1; (pi — d\ — 1)N + 2 + di] of Eq. (l37I) . namely 

71 )6l _ dl (A)71 t2 +di(/^) = — 1 — L ^ )2 +d 1 (/i)7"i ibl _ ( i 1 (A) 

+ E Lf , + uT T 1 ,g(/i)Ti, 2+bl _ g (A) ) for6i>iV. 

e=2+ bl -N 



(38) 



On the other hand for b\ < N we have to make a combination of two entries [1; (&i — d\ — 
1)N + 2 + di) and [1; biN + 1] of Eq. fl371) . These equations can be combined to produce a 2 x 2 
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Operator 


Creation 
index 


Combination 
index 


Number of 
equations 


%.,bx-di (A)^l,2+di (aO 

ai = 2 


bi > N 


c = b\ — d\ — 2 


1 


^l,6i-di(A)^l,2+di(/^) 

a x = 2 


bi < N 


c = b\ — d 1 — 2, bi — 1 


2 


TlM-di {fl)%,i-l,ai+di (A) 

3 < ai < iV + 1 - bi 


h<N 


c = 0, ...,6i - 1 


&1 


%.,h.-di \l*)1ai-l,a\+di (A) 

iV + 2-6i<a 1 </V 


bi < N 


c = 0, . . . , TV — ai 


JV + 1 - ai 


%.,bi-di {^)%ix-l,at+di (A) 

3 < ai < N 


h>N 


c = bi — N, . . . ,N — ai 


2N + 1 - oi - h 



Table 2: The linear combination derived from Eq. (j37|) to achieve appropriate suitable commutation rules between 
the fields 7i,b 1 _ ( j 1 (A) and T ai ^i^ ai+c i 1 W- The indices a\, b\ and d\ belong to the intervals defined by Eq. 



system of linear equations given by, 



A 



(budi) 



(A,/i) 



71,1(^)71^+1 (A) 
71,2(^)71,6! (A) 



i2(A, //)};} 



%.,bi-di (A)7i,2+d 1 (/i) 

71,1+6! (A)Ti,i(//) 



6i+l 



\bi-di,2+di 



~ £ ( ^'^f+iT I ^ri^CA) for 6i < TV, 
J R(A,/i) 2 1 +bl Lg i g 



e=3 
(ai,6i). 



where the 2x2 matrix A 3 1 ' (A,//) is 



.4 



(6i,di) 



(A, A*) 



/«A./ ( );;::':r"' : />'( a. // )!' : 



'6i,2 
\6i+l,l 



(39) 



(40) 



^(Aj^bi+lj -R^aO^ 1 

The desirable commutation relations are derived by solving Eq. (l39l for the product 7^ 2 (a0 
7i,6i(A) with the assistance of Cramer's rule. As a result we obtain a single linear equation 
relating the terms 7i j 6 1 _d!(A)7i i2 +di(A*) and Ti : 2+d 1 (p)%.,bi-di(^) B whose solution gives us the 



We recall that the product / 71,2+di(A t )^i,6i-rfi(A) appears in the sum of Eq. ([551) for e = 2 + di < b\. 
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following expression, 



-di W^l,2+d 1 (aO 




b\-d\,2+d\ 
b 1 -d 1 ,2+di 

6l+l,l 



b l -d 1 ,2+d 1 



^1,2+^1(^)^1,61-^1 



(A) 




R(\, fi){jR(\, H) b b ] 



+ 



fei+i 



E 



R(X,fi) 




R(X,fi) 



b 1 -di,2+d 1 
2+6i -g,e 

72+6i-e,e 



^l,e(/^)^l,2+6i-e(A), 



e=2 



J R(A,/i) 1 ; 1J R(A,/i) 




e^2+di 



for 61 < AT. 



(41) 



Let us now discuss the strategy we have used so far in order to obtain the commutation relations 
for the creation fields. The basic idea is that any possible product of creation operators on the 
right-hand side of the commutation relations should be equally ordered as far as the rapidities A 
and fi are concerned. This ordering is certainly the opposite of that we have started with for the 
left-hand side product of creation operators. In addition to that, the specific ordering choice is 
actually dictated by the results established in section 13.11 for the commutation rules among the 
diagonal fields and the basis vectors. The commutation relations for the creation fields are then 
constructed to bring a given creation operator with the transfer matrix spectral parameter to the 
further left position in products of monodromy matrix elements that are not proportional to the 
eigenvectors. 

The above procedure is of special importance when we deal with the commutation rules be- 
tween the operators Ti bl _ dl (n) and Zn-i,ai+di(A) for 3 < a\ < N. In this situation the right- 
hand side of the respective commutation relations is not allowed to possess terms of the form 
/ 7i,6i-e(A i )'7ai-i,ai+e(A) with e > 0. This can be clearly seen by considering in an explicit way the 
linear systems associated to such commutation rules. From Table 2 we observe they depend very 
much whether bi < N or b\ > N. The former case is rather similar to the linear combination 



17 



discussed in previous subsection and we find that 

-l,ai+6i— 1 (A) \ 
-l,ai+6i— 2 (A) 



4 ai ' 6l) (A,M) 



ai— 1 



e=l 



\ ^l,6i(^)^i-l,ai(A) / 



^g,ai+l(A)'^ ll -e,fei^l(A i ) 

y ^i,a 1 +6 1 -i(A)'^ ll _e,i(A i ) y 



oi+fti— 1 



E ^ ai ' bl) (A,/i)T 1>e -(/i)T ai _ liai+bl _ e -(A) for6 1 <iV and 3 < fll < iV + 1 - 6j (42) 



c=6i + l 



4 01,6i) (a,m) 



^l,oi+bi-iv(A t )^ai-l,Ar(A) 
^l,ai+6i-iV+l(A t ) / ?ai-l,Af-l(A) 



V 

/ % jai {\)T ai _ EM {^) \ 



ai— 1 



E^/C-S 



e=l 



X 



^1,61(^)^1-1,01 (A) y 

i-e,6i(/^) 
^g,ai+l(A)'7^ 1 _e i f, 1 -l(/i) 



A? 



4 ai,6l) (A, ) u)Ti,e(//)7; i _i jai+6l _g(A) 

e=6i+l 



I ^g,Af(A)Zii-e,ai+bi-Af(A t ) J 

for bi < N and JV + 2 - ^ < a x < iV 



(43) 



where the b\ x &i matrix A^f 1,bl \\, //), the 61 x 1 vector v{ ai ' bl \\,n), the (JV + l-ai) x (iV + l-a x ) 
matrix A^ ai ' 6l) (A, /x) and the (Af + l-ajxl vector u$ ai ' 6l) (A,^) are given by Eqs. (I22lf23lf25lf26|) . 
respectively. 

For b± > N we have instead a different linear system of equations. In this case the relations 
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coming from Eq. (137|) for c = b\ — N, ■ ■ ■ ,N — a\ are organized as follows, 

^l,ai+6i-7v(A t )^ 1 -l,Ar(A) ^ 
^1,01+61-^+1(^)^1-1,^-1 (A) 



4 ai ' M (A,/i) 



ai — 1 



e=l 



X 



Y ^l,7v(/^)Zii-l,ai+6i-Af(A) y 

,ai+6i-JV 
,01+61 -JV+l (A)^j 1 _g i AT_i (/i) 



for bx> N and 3 < a x < iV, 



(44) 



\ ^g,iv(A)^ai-e,ai+6i-jv(/ i ) y 



where the (2N + 1 — ai — 61) x (2N + 1 — ai — 61) matrix A± 1 (A, /i) is given by 

p/\ \ai+6j-Af,AT R / X aoi+6i-AT,7V r>(\ \a 1 +b 1 -N,N 

^K^i H')N,a 1 +b 1 -N rL V A ) WaT-1,oi+6i-JV+1 "' ^H'v /^ai+fii-i^iV 



.4 



(ai,6i) 



(A,/i) 



D / \ \ai+6i-Af+l,AT-l p /\ „\oi+6i-Ar+l,Ar-l n/\ no 1 +& 1 -AT+1,JV-1 

"V^j t Jj )N,a 1 +b 1 -N rl V A ) r^iV-l,ai+6i-iV+l '" "VA h i )a 1 +b- l -N,N 



U{\ ,.\N,ai+bi-N td(\ ,,yv,ai+oi-jv p/\ yv,ai+o ; 

\ - K l A )A t JjV,a 1 +bi-AT K V A >/ i JjV-l,ai+6i-iV+l '" " K l A ' A t J ai+fol _ 



^,01+61-^ 



^7V,ai+6i-JV 
-N,N 



J 



(45) 

Direct inspection of Eqs. (I42II43||44|) reveal us that the variables of the corresponding system of 
linear equations are indeed the products / 7i jfel „ dl (/i)7^ 1 _ l ai+dl (A). The origin of these independent 
linear systems is directly related to the existence of three distinct intervals for the index d\ once the 
variables a\ and &i are fixed. Their structure have been constructed to collect together products 
of creation fields whose respective values for the index di belong to one of the three possible such 
intervals. The commutation relations between the operators Ti^-^in) and T^-i^+^A) are 
then determined by means of the systematic application of Cramer's rule to solve the systems of 
linear equations (I42|43II44I) . The procedure is similar to that already described in section (13. ip . 
For instance, we note that the reversed product terms appear on the first non-homogeneous part 
of Eqs. ( l42ir43lT44l) when the sum index is e = a± — 1. The task is however much more cumbersome 
since we have to apply Cramer's rule for each left-hand side product of creation fields entering 
Eqs. (142(431441) . The number of different commutation rules associated to each linear system for 
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a given multiparticle state depends strongly on the index b\. This feature and the dependence of 
the total number of relations for a given N is illustrated in Table 3. 



Linear 


Number of Commutation Rules 


fixed bi 


fixed N 




AT+l-6i 
ai=3 

= (N — 61 — 1) (&i — 1) 


X^(N - bi - l)(bi - 1) 

bl (N — 1)(N — 2)(N — 3) 
fi 




N 

E ( N + 1 -^) 

oi=AT+2-6i 

_ (6l " l)6l 
2 


v^ 1 (h - l)h 

^ 2 

61=2 

N(N — 1)(N — 2) 
fi 




2JV-&1 

E (2iV + 1 - ai - 61) 
(2N -b x - l)(2JV-6i-2) 


2 ^ 2 (2iV - b x - 1)(2N - 61 - 2) 


^ 2 

bi=N 

N(N — 1)(N — 2) 


2 


fi 



Table 3: The number of commutations relations between the creation fields for fixed &i as well as 
for fixed N. 

We would like to conclude by presenting an example of a complete set of commutation rules 
between the creation operators. The simplest non-trivial situation occurs when the respective 
azimuthal spin component is indexed by b\ = 2. The intervals (1361) tell us that such relations are 
those between the operators T lt2 (fi) and 7^ 1 _i >01 (A) for 2 < a x < N. The case ai = 2 follows from 
the general expression fHTj) while for the remaining values we need to solve the linear systems 
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( 14211431) . Putting the results together we find that, 




ai,e-2) 



0i,2 

oi+l,l ig . a i 

01+1,1 



(A, ^)7i 1 g(/i)^ 1 -l,ai+2-e(A), 



l(A)T a 



oi+l,l 
ai+1,1 



ai,2 
fti+1,1 



for 2 < a! < - 1 



(46) 



^l,2(/i)^7V-l,Af(A) 



-R(A, /^)^2 




^g,iv(A)7^v_e,2(Ai) 




^l,e(/ i )7/V-l,AT+2-e(A), 



(47) 



where we recall that function (A, //) is given by Eq. (l30|) . 

3.3 The creation and annihilation fields 

The third class of commutation rules that still need to be considered are those between the basis 
vectors ^i,b(/x) and all possible annihilation operator T a ^+d x ,ai-i(A). We stress here that the indices 
ai, d\ and h are the same used in subsection (13. 2p and therefore they belong to the intervals defined 
by Eq. (l36|) . In general, the result of taking the operator T ai+di ai _i(X) through the creation fields 
Ti^{fi) is the generation of several distinct annihilation fields besides T ai+ d ltai -i(^)- The basic 
strategy we shall use to construct these relations is to enforce that all such resulting annihilation 
fields must be kept in the right-hand side position on the respective commutation rule. A detailed 
study of the Yang-Baxter algebra reveals us that these commutation rules can be built up from the 
entries [ax+di+ci(N— 1); (g^ — 2)A^ + 6+c 2 (A^ — 1)] of Eq. ffTTl) . The structure of the corresponding 
equations are now more complicated since they carry a dependence on two independent indices C\ 
and C2 as well as on the particular combination a\ + d\ . Taking into account the latter feature we 
find convenient here to express the expressions by means of the following auxiliary indices, 



fi = a i + di and f 2 = b — di — 2 



(48) 
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The starting point relations are then obtained from Eq. (TTTl) by choosing a = ci + 1, b = fx — cx, 
c = ax + C2 — 1 and d = b — C2, namely 
h 

2j R(\, ^)//i^ Cl +i / 7g,ai+C2-l(A)7} 1 +l-e,6-c 2 (A t ) = 
e=l 

m{ai+b-2,jV} 

£ r cl+lie -( / x)r /l _ cliai+6 _ 1 _ e -(A)i?(A, / x)^+ri^-7 2 - ( 49 ) 

e=M{l,oi+b-JV-l} 

The analysis of Eq. (|49|) turns out to be far more involved than that discussed in the two last 
subsections. The respective linear combinations have to be performed considering two distinct 
steps and they will culminated in two coupled linear systems. It is rather illuminating to start 
this discussion by examining Eq.( l49l) in the particular situation c\ = C2 = 0, 

h 

^j2(A,Ai)S5 +1 - e r 8l «- 1 (A)r A+1 _ a , 6 ( M ) = 

e=l 

m{ai+fe-2,V} 

7l >g (//)T fliai+b _ 1 _e(A)ii:(A,//)^;;i b 1 _e,e> ( 50 ) 

e=Af{l,oi+6-JV-l} 

Let us now concentrate our attention on the left-hand side of Eq. fl50l) . For e = fx it gives us the 
product 7^ 1+ d ll a 1 -i(A)7i i 6(/i) which is indispensable to build up the desirable commutation rules. 
Unfortunately, for arbitrary values of b, this term is not the only one involving annihilation and 
creation operators generated by Eq. (l50l) . From the left-hand side of Eq. (l50l) we see indeed that 
7^ ai _i(A) behaves as annihilation or diagonal operator when e> a\ — l while T ai+dl+1 _ sb (fi) play 
the role of creation fields for e > ax — fi = Oi — (b — dx — 2). This product combination clearly spoil 
the main characteristic one would expect from the commutation rule among 7^4.^,01-1 (A) and 
Tx^(fi). The first task therefore is to eliminate the left-hand side products 7^ iai _ 1 (A)7^ 11+(il+1 _g i b(/^) 
when the index e takes values on the intersection of the above intervals. Considering the upper 
value of e in Eq. ([5D|) it is not difficult to see that the mentioned interval is 

a x - m{l, f 2 = b — dx — 2} < e < fx = ax + d x . (51) 

The constraint fl5T|) tells us that the case b = 2 is a fortunate exception and therefore there 
is no need of considering the above mentioned cancellations. This means that the appropriate 
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commutation rules between T ai+ d 1)ai -i(ty and Tl^O") can be derived directly from Eq. (l50l) . After 
few manipulations we find, 



7X,ai-i(A)7i )2 (^) 



Ti j2 (//)7> 1)ai _i(A) + 2^ ^ t r \ s ,\/i,i e ' e r i^(^)^/ 1 , ai +i-g(A) 



01 



-E 



— r^^CA)^^^). (52) 



On the other hand, for b > 3 we are forced to get rid of the products 7^ jai _ 1 (A)7^ 1+dl+1 _g i b(/x) 
when the index e belongs to the interval (jBTj) . This can be done by exploring the linear combination 
of equations derived from Eq. (j4*9j) for suitable values of the index c\ while c 2 is fixed at c 2 = 0. 



The structure of the linear combination depends on the sign of the azimuthal spin of the product 



x ox-i(A)^i as ^ is illustrated in Table 4. 



Operator 
eliminated 


Creation 
index 


Combination 
index c\ 


Number of 
Equations 


'Tg.ai -1 (A) Tft +l-g,b ( aO 


b-2<d x 


ci = 0,...,6-2 


6-1 


^g,ai-l(A)'7} 1 +l-g,fe(A i ) 


b-2>d l 


ci = 0, . . . , c?i + 1 


rfi + 2 



Table 4: The linear combination from Eq. (|49|) for C2 = to cancel the operators 75 iai _i(A)7/ 1 +i_e,b(/u) when e 
belongs to the interval (|5Tj) 



Taking into account Table 4 it is not difficult to start the construction of the first needed linear 



Recall that from Eq. (fT5| the corresponding spin value is /a = b — 2 — d\ 
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combination for b > 3. Its main structure is as follows, 



and 



^ ^/i-b+2,ai-l(A)T 6 _ 1;6 (/x) ^ 
^/l-6+3,ai-l(A)^_ 2 , 6 (/x) 



ai+6-2 



e=l 



X 



( Ti ig (/x)7) 1)ai+6 _i_e(A) ^ 

^2,e(^)3/i-l,ai+6-l-e(A) 



\ ^6-l,e(/^)7fi-6+2,ai+6-l-e(A) y 



/i-6+l 

E 4 /l ' fe) (A, /U )T g , ai _ 1 (A)T /l+1 _ e - b ( /U ) 

e=l 



for 6-2 <di, 



4 ai,dl) (A,/i) 



/ 



x 



^ ^i-l,ai-l(A)7^ 1+2 ,b(/^) 
^l,ai-l(A)7d 1+ l,&(/^) 



\ ^>i,ai-l(A)71 i6 (^) ) 

\ 



m{ai+b-2,N} 

E ^/C;^- 

e=M{l,ai+6-A r -l} 



^l,e(A t ) 7 /i,a 1 +6-l-e(A) 



^2,e(/^)2/l-l,ai+&-l-e(A) 
\ % 1 +2,e(^)%i 1 -l,a 1 +b-l-eW J 



ai-2 



^ E 4 ai ' dl) (A,/i)T g , ai _ 1 (A)T /l+1 _ e -, b ( / u) 



e=l 



for b-2>d 1 . 



(53) 



(54) 



The (6 - 1) x (b - 1) matrix A { 5 fub) (X, /i) as well as the (^ + 2) x + 2) matrix ^° 1,ai; (A, a«) 



(oi,di) , 



are given by 



RiX^Y/J^- 1 R(X, fi) j f 1 ~ 1 b+3 ' b ~ 2 



-6+2,6-1 -6+2,6-1 / 



/i -6+3,6-2 



(55) 
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and 



A 



(oi,di) 



(A,//) 



R{\^)lf +1 



R(\») f fli 
R(\,rijlt, 



(56) 



while the (b — 1) x 1 vector v^ 1,b \\,fi) and the (di + 2) x 1 vector VQ ai ' dl \\, fi) are 



(A,//) 



Ri^vYff-ii 



\ 



\R(X,fi)% f itl2,li J 



( R(\nr f (\ +1 - B ^ 



(ai,dl) 



(A,//) 



e,/i+l-e 



V^(A,^ti;i 2 ; 



(57) 



By construction the linear systems defined by Eqs.( l53|54l) contain explicitly the product 
/ 7j i ai _ 1 (A)7^b(/i) which can be calculated by means of Cramer's rule. From the right-hand 
side of Eqs.f l53|l54l) we observe that this solution is able to produce products having the gen- 
eral form 7^ 1+ i j g(w) / 7/ 1 _ Cl , ai +b-i-e(A). This type of terms include the desirable reversed product 
T\fr{\L)Tf x ,ai-i(A)|j but also several other products having rather undesirable properties. The lat- 
ter feature occurs when e < c\ + fi = b + c\ — d\ — 2 because the operator 7/ 1 _ Cli(ll+ ;,_i_g(A) 
turns out to be a creation field and its presence on the further right position of the right-hand 
side commutation rule should definitely not be permitted. We stress here that this problem is 
independent of the character of the companion field 7^ 1+ i i g(/i). This is clear when 7^ 1+ i i g(/i) plays 
the role of either annihilation or diagonal operator since we have exactly the same situation we 
managed to handle for the left-hand side of the commutation rules. The case when 7^ 1+1 g( / u) 
play the role of creation operators is more subtle because the product 7^ 1+1 g(/i)7/- 1 _ ci ai+fo „ 1 _g(A) 
could, in principle, contribute to the eigenvector basis. The order of the rapidities in such product 
is however opposite to that already chosen in section (13.21) . For this reason such remaining type of 
product combination has also to be avoided. Therefore, no matter the role of 7^ 1+lj g(/z), a second 
step is still necessary to provide us the means to compute the product 7^ 1+ i i g(/i)7/ 1 _ Cljai+ ; ) _i_g(A) 



4 This is the first component of the right-hand side vector (c\ = 0) in Eqs. (l53|54p when e = b. 
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as long as the index e satisfies the relation, 

e < ci + f 2 = b + ci — di — 2. 



(5f 



This additional task has to be implemented without spoiling the main construction underlying 
the first linear system. This can be done by using the freedom of the extra index c 2 since so far we 
have kept it fixed at c 2 = 0. In Table 5 we have summarized the second linear combination for a 
given C\ with the help of index c 2 . The equations derived from such linear combination are suitable 
to calculate the products 7^ 1+lj g(/i)7} 1 _ cliai+ { ) _ 1 _e(A) on the interval M{1, a\ + b — N — 1} < e < 



m{ai + 6 — 2, N}. Recall here that this interval is the total range of the index e in Eqs.( j53ll54l) 



Operators 


Creation 


Combination 


Number of 


calculated 


index 


indexes c 2 


equations 


^i+l,e(^)^/i-ci,oi+6-l-e(A) 


b-2<N-a 1 


c 2 = di — ci + 2, . . . , b — 1 


ci + b - d x - 2 


%i+l,eWTf 1 - Cl ,a 1 +b-l-e{X) 


b - 2 > N - ai 


c 2 = d\ — ci + 2, . . . , 


Ci + N — ai — di 






N - ax + 1 





Table 5: The linear combination from Eq. (|49|) to compute products of the form 7^. 1+ i : g(/i)7} 1 _ Cl . ai+ b_i_g(A) 



By substituting the data of Table 5 in Eq. fl49|) we find that the second system of linear equations 
are given by, 



(/l-Cl,/ 2 +Cl) 



(A,/x) 



Zi 1 +l,l(M) 7 /i-ci,a 1 +6-2(A) ^ 
^1+1,2(^)^/1 -ci,ai+6-3(A) 



e 

+1 



e=l 



\ ^l+l,/2+Cl(^)7/l-Cl,/l-Ci+l(A) J 

^ ^g,/i-ci+i(A)T fl+1 _g i/2+Cl (/i) ^ 
/ ?i,/ 1 -c 1 +2(A)T /l+ i_g i/2+Cl _i(^) 



ai+6-2 



v. 



(/l-Cl,/2+Ci) 
7 



(A, A*) 



e=fe+ci — di—1 



y ^g,oi+6-2(A)7/ 1+ i_e i i(/i) y 

1— ci,ai+6— 1— 

.(A) 



for b < N - a x + 2, 



(59) 



2(3 



and 



|(/l-Cl,/2+Cl) 



(A,/i) 



^ ^ 1+ i j0l+6 _jv-l(/i)7fi-ci,iv(A) ^ 
^i+l,ai+6-Jv(/^)7fi-ci,iV-l(A) 



X 



\ ^l + l,/2+Cl(/^)7fl-Cl,/l-Cl + l(A) J 

^ ^,/i-ci+l(A)7fi+l-e,/ 2 +ci(^) ^ 
/ 7g,/l-ci+2(A)7} 1+ i_g i / 2+Cl _l(/x) 



— e 
,ci+l 



e=l 



\ 7i,jv(A)7} 1+ i_g iai+fe _Ar„i(^) y 
X / ^ 1 +l,e(A t )7/i-ci,ai+6-l-e(A) 



- E 

e=b+ci— di — 1 



for 6 > TV - ai + 2, (60) 



where the 6x6 matrix ^ a '^(A, yu) and the N + 1 — a x N + 1 — a matrix A^' b \X, /i) are given in 
terms of the structure constants as 



4^(A,/i) = 



\a+b,l 



a+b,l 



(61) 



and 



4 a ' b \\^) 



N,a+b-N+l N-l,a+b-N+2 



1,6 



-1,6 \ 



R(\ A0jV,a+6-iV+l R(\ A*) N-l,a+b-N+2 



a+2,fe-l 



R( X i ^)a+l,b 
R(\v) a atli 1 



I D/\ n \N,a+b-N+l R / X A N,a+b-N+l p/ \ , \iV,a+6-iV+l 

y -K(,A,/ij Ara+{ ,_ Ar+1 n{A, fj,) N _ la+b _ N+2 ■■■ n{A,fi) a+lb 



>N,a+b-N+l 



(62) 



while the 6x1 vector ^'^(A, ji) and the iV+l-ax 1 vector ^'^(A, /i) are 



(a,6). 



( Til \ „^+l,i 



4 a ' b) (A,/i) = 



R{\v) a a 



+6+l-e,e 

a+2,6-1 
a+6+1— e,e 



( r>i \ i<H I.. 6 



(A,p) = 



i?(A,/i) 
i?(A,/x) 



a+6+1— e,e 

a+2,6-1 
a+6+1— e,e 



\ ^lA,A t Ja+6+l-e,e / 



(63) 
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Annihilation 
index 


Combination 
index 1 


Combination 
index 2 


b - 2 < d x 


ci = 0,...,6-2 


c 2 = dx — cx + 2, . . . , b — 1 


b-2> d x b - 1 < N - ai 


cx — 0, . . . , dx + 1 


c 2 = dx — Cx + 2, . . . , b — 1 


b-2> dx b-l> N -ax 


cx — 0, . . . , dx + 1 


C2 = c?i — ci + 2, . . . , N — a x + 1 



Table 6: The linear combination from equation P5| necessary to obtain suitable commutation rules between the 
fields T ai+duai ^i(\) and T ltb (fj,) 



Let us now show how the first and the second systems of linear equations work in practice 
together. In order to see that we shall discuss in detail the commutation rule among the fields 
^/i,ai-i(A) and 7i i3 (/x). We start by considering the first linear system of equations (E3E1D which 
for 6 = 3 gives us the following expressions, 

,(/,.:■!),, 1 ^/L-l,oi-l(A)^s(^) \ , sai-1,3 / ^l,e(/^)^/i,ai+2-e(A) 

^5 (KM I | = 2^ i? ( A ^)ai+2-e-g 



T fljai _x(X)T lj3 (fM) 
/l-2 



e=l 



^2,e(A t )^/i-l,a 1 +2-e(A) 



-E 

e=l 



,(/i.3) 



(A, /i)Tg iai _ 1 (A)r /l+1 _g >3 (^) 



for tZi > 1, (64) 



and 



4 G1,0 W) 



Xn-i,ai~i(A)3^ 3 (/i) 



^l,e(A*)^oi,oi+2-e(A) 



m{ai+l,AT} 

e=M{l,ai+2-AT} \ ^2,e(A f )'?a 1 -l,ai+2-e(A) 



7"a 1 ,a 1 -l(A)T li3 (/i) 

ai-2 

: ^(A^)ai+2~ e -e - ^ (A,/^)^ ! a 1 -l(A)7; i+ i_e,3(/^) for = 0. 



(65) 



e=l 

The matrices A^ 1 (A, /i) and A 6 a (A, /x) are given by 



4 /l,3 W) 



/?(A,/i)^:g ^(A,//)^ 



4 G1,0 W) 



-R(A,/i)^_i;2 -R(A, ^)oi-i,2 



while the vectors i^'^A, /x) and Ug ai,0 ^(A, //) are 



(66) 



(67) 
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By using Cramer's rule in Eqs. (I64f65l) we can compute the products 7/ i ai _ 1 (A)7 Li 3(/i) and the 
results are, 

ai+1 a(a, ri:\zil e , e R(\, ^)fc 1,2 



/i,ai-l 



E 

e=l 



oi+l 



- E 



e=l 



/i-2 

E 

e=l 



r(\, fi) f ;-{' 2 2 r(x, nY f it lj2 



l,ai+2-e 



(A) 



R{\li)%-^ R(X,^ 



fiA 



r(\^ 2 r(\, »y£ +1 -~ e 



ai+2-e 



(A) 



R{\^y f \Z\l R(\,v) e fX 



e,/i+l-e 
2 



^(A,/i)fc 1,2 



i2(A,Ai)g, a 
Z(A,//)£ 
for di > 1 



fl(A, J2(A,/i)£i, a 



^g,ai-l(A)7 fl+ i_e,3(/^) 



(6* 



and 



(A)T li3 (/i) 



m{ai+l,Af} 

E 

e=M{l,a 1 +2-A r } 



m{ax+l,iV} 

E 

e=M{l,ai+2-AT} 



ax -2 

E 

e=l 



-R(A,/i)^ + 2^g ig -R(A,/i) 



^(A,//)^_^2 -^(A, ^) ax-1,2 



^(A,/i) 



ai-1,3 p/-v nox— 1,2 
ai+2-e,e- rl l A ! Pjaj-1,2 



i2(A,/i) 



'ax,l 

ax-1,2 
ax-1,2 



R(\^r i 



R(\,li) 



CL\ ,1 

ax,l 
ax-1,2 



^A,A*)S^ a i2(A,/i)S;i +1 " e 

^(A,^)oJ-li2 ^(Aj^af-^ 



ox — l,ai+2— e 



(A) 



ax, ax +2— e 



(A) 



'IX, 1 



i?(A,/i) ai1 



^(A,/i)aJ_iJ -R(A,//)aJ-i,2 



ax,l 
ax ,1 



% 



e,a\ — 1 



(A)T B 



ax+l— e,3 



for = 0. 



At this point we note that the commutation rules 



between 7/i=ai+di,ai-i 



(69) 

(A) and 71 i3 (//) 



possess the appropriate form only when d\ > 1. In these cases the right-hand side of Eqs. 
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does not generate undesirable products of fields and the second system of linear equations is 
not needed. This is not the situation of the remaining values d± = 1 or d% — 0. For example, 
when d\ — 1 we see that the first term of Eq.( l68l) for e = 1 produces the undesirable product 
T 2j i(n)T aijai+ i(X) . This product can, however, be eliminated thanks to the second linear system of 
equations. Indeed, by considering Eq. fl59|) for c x = 1 we can compute the product 7^ 1 (/i)7^ i ai+1 (A) 
as, 



7i,l(/^i,ai+l(A) 



ai+1 R(X /i) e ' ai+2_e 

ai+1 Rl'X //"\ ai+1,1 



pT 2ie -(/i)T aiiai+2 _ e -(A) for a x < N - 1. (70) 



g=2 ^(^'^)ol+l!l 

For the value d\ = 1 we can now substitute Eq. flTOj) in Eq. fl68|) . This cancels out the undesirable 
product, providing us the suitable commutation rule between the operators / 7^ 1+lai _ 1 (A) and 



- i ai+l,ai-l 



{mM = - E 



+1,1 



ai+1 j-,/ \ ..\oi,2 
, _2 ^(A; AfC+1^1 



-R(A,aO"1 + 1;'i 



^(A,^)aJ+l',l 



#(a,aC 



ai + 1,1 



+2-e,e 



^(a,ac 2 



+1,1 

\ai,2 



ai+1,1 
ai,2 



T 2,e(^) 



X Zn,ai+2-e(A) 



ai+1 

E 

e=l 



ai+2-e,e oi,2 



oi,2 



^(A, A*)aJ+l,l ^(A,//)ai^l 



^l,e(/^)Zii+l,ai+2-e(A) 



Oi+l 



+ E 



e,ai+2— e 



^aC+V^aC 



ai— 1,3 
+1,1 



e=l 



^(A,aC 



Ol+l,l 

+1,1 



^ai,2 
u+l 

lOi ,2 



ai— 1 



- E 



^(A,A*)aJ+l,l 



e,ai+2— e 
ai+1,1 

e,ai+2— e 



^g,Oi+l (A)^!+2-e,l (/^) 



#(A, 6 



e=l 



-R(A, A*)oi+i,i ^(Aj A*)ai+i^i 



ai+1,1 



ai+1,1 



^i,ai-l(A)^j 



ai+2-e,3 



(A*) 



(71) 
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The case d% = is more laborious due to the presence of three different types of undesirable 
products. They come from the first (e = 1,2) and the second (e = 1) terms of Eq. (!69|) and 
their respective forms are T 2 ,i(ii)T ai - hai+ i(X), T 2j 2(fi)% 1 -i, ai (X) and T 1A (fi)T auai+1 (X) . All these 
terms can be eliminated with the help of the second linear system of equations as follows. The 
last undesirable term T L1 (/i)7^ i ai+1 (A) is the simplest one to be computed. Its expression comes 
directly from Eq.( |59l) by taking the values c\ = d\ = 0, 

7l,l(^)Zii,ai+l(A) = ^ — — ?ai+ll ^,oi+l(^)^i+l-e,l(^) 

g=1 K {W) ai+ i,i 

ai+1 d/\ \oi+l,l 



fl(A,A*)« 



e=2 



S|=fr 1)e -( M )r 01)ai+2 _ e -(A) (72) 



The remaining products 72 i i(//)7^ 1 _i )01 4.i(A) and 72,2(^)^*1-1,01 (A) are obtained as a solution 
of a 2 x 2 linear system of equations. It is derived from Eqs. fl59ll60p for the values C\ — 1 and 
d\ = and the final result is, 

,(ai-l,2) M x| ' ? 2,l(Ai)7; i _i >ai+ i(A) \ g a +1 _g / 7i iai (A)7; i+ i_e,2(/i) 

\ 72,2(^)^-1,01 (A) / e=l \ 79 |01+1 (A)7^ 1+ i_e,i(At) 

Ol+l 

-I]4 ai " 1,2) (A,/i)72, e -(/i)T ai _ liai+2 _ e -(A) for a x < N - 1, (73) 

e=3 



A? 



-R(A, /i) Ar | 2 7 2)2 (^)77V-l,Ar(A) = -R(A, /i)^_ 12 e 7g i Ar(A)T/v + i_g l 2(A t ) 

e=l 

AT 

- ^ J R(A,/i)^ 2 _g i gT 2i e(/i)T 7 v_i,iv+2-g(A), (74) 



e=3 



where 

4^-1,2)^ )= *(a,aos3 ] and V ( 01 -1,2) (A )= ^(A^)::v 2 2- e -, e - 

V^(A,^S1 i2(A,/i)SS 1 ' 1 / ' \«(A, 

(75) 

By applying Cramer's rule in Eq.( |73l) we can easily compute the undesirable terms 
72,i(//)7^ 1 _i )0l +i(A) and 7 2i 2(/w)7^ 1 _i jai (A) for a\ < N — 1. We then substitute this result as well 
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as the expression (!72|) for 7i i i( / u)7^ ljai+ i(A) in Eq. (!69|) . In this way we are able to obtain the 
commutation rule for the fields 7^ liCll _i(A) and 7i )3 (/^) for a\ < N. For a\ = N we just have 
to eliminate the term ^^(/^T/v-i.aKA) of Eq. (l69l) with the assistance of the relation (17^1) . The 
final step to obtain suitable commutation relations concerns with the reordering of the products 
T aii£tl (A)T lj2 (» with the help of Eqs. fl28ll29|) . 



4 Identities among the weights 



In this section we describe a procedure to derive identities between the amplitudes -R(A, /i)^ from 
both the unitarity relation and the Yang-Baxter equation. These identities are essential to carry 
out simplifications on the eigenvalue problem without the necessity of referring to specific weights. 
We shall start by considering the consequences of the unitarity relation ([5]). 

4.1 Unitarity relation 

A systematic study of the unitary relation ([5]) revealed us the existence of two key sets of inde- 
pendent weights. They will play a relevant role in the analysis of Eq. ([5]) and are defined as, 

j>qi _ I R{\tig£% for j = l 



alf = { — (76) 

\N-q 1 +c,N+l-c 
>N+l-b,N-qi+b' 



Kf = { ( 77 ) 

R(^ N+V-bJf-qi+bi for J = 2 ' 

where q% = 1, . . . , N and b, c = 1, . . . , q\. 

The above definition explores the block form of the i?-matrix on the basis of the U (1) operator 
S z ® In + In ® S z . The index j = 1, 2 is used to split the charges with q = 1, . . . , N (j — 1) 
from those with q = N + 1, . . . , 2N — 1 (j = 2). Considering Eqs. fl76|77l) we are able to write the 
matrix elements of unitarity relation (jSJ) as, 

U[b, c]f = <S4 9 c - $b,c = for qi = 1, . . . , JV; b, c = 1, . . . , g l5 j = 1,2. (78) 

k=l 
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A non-trivial identity between the weights is already relevant for the analysis of the one-particle 
problem. It is obtained from the 2]f component and by substituting the respective indices in 
Eq.(EED we find, 

The type of relations are necessary because the combined use of the commutation rules, dis- 
cussed in previous section, forces us to reorder the spectral parameters of some i?-matrix elements. 
It is fortunate that we can solve this problem with the assistance of the unitary property. In gen- 
eral, more complex identities are needed and they are derived by analyzing linear systems of 
equations based on Eq.( |78l) . The first such system is obtained by means of the linear combination 
U[b; i + l]} +2 a.i'*j+2 ~ ^[^> * + 2]} +2 ffi'^ for b = 1, . . . , i + 1. These equations make it possible to 



write the following linear system for variables a 



-j,i+2 -j,i+2 -j,i+2 -j,i+2 



b+l,i+l a 'l,i+2 



a b+l,i+2 a l,i+H 



( j,i+2 j,i+2 



1,2 



*2,2 



j,i+2 
\ a i+l,2 



1,3 



J,i+2 i,i+2 



l 2,3 



a: 



j,i+2 
i+1,3 



i,i+2 \ 
1,i+2 

J4+2 
^2,1+2 



a i+l,i+2 / 



( -j,i+2-j,i+2 -j,i+2-j,i+2 \ 



a 2,i+l a 'l,i+2 



a 2,i+2 a l,i+l 



-j,i+2 -j,i+2 -j,i+2 -j,i+2 



\ 



a 3,i+l a 'l,i+2 



-j,i+2 -j,i+2 
a i+2,i+l a l,i+2 



Q 3,i+2 a l,i+l 



-j,i+2 -j,i+2 
a i+2,i+2 a l,i+l J 



( 



-j,i+2 
\ a ti+2 ) 



The solution of this system of equations for three particular components is able to generate 
classes of important identities. The first family is found by solving the system (IHUjl for the first 
component. By employing Cramer's rule we found, 



j,i+2 
<2 



J,i+2 
l i+l,2 



a: 



j,i+2 
Z l,i+2 



j,i+2 
'i+1,1+2 



d,i+2 -j,i+2 



1,i+l 

-j,i+2 
Z 2,i+1 



1,1+2 

-J,i+2 
l 2,i+2 



a l,i+2 ( 1 1 



j,i+2 

h,3 



j,i+2 
<3 



J,i+2 
l l,i+2 



j,i+2 
l i,i+2 



The remaining families of identities follow from the solution of Eq. (l80|) for the second and last 
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components. After making the ratio of these solutions with Eq. llHTT) we have, 









j,i+2 j,i+2 
a l,2 °1,4 


■ ■ ■ 


+2 
+2 


-j,i+2 
a l,i+l 


-j,i+2 
a l,i+2 












-i,i+2 
a 3,i+l 


a 3,i+2 




j,i+2 j,i+2 
a i,2 a i,A 


Q j,j+2 


-i,i+2 
a l,i+l 


a l,i+2 






i,«+2 

«1,3 • • • 


i,«+2 




a j ' i+2 


a 2,i+2 






7,*+2 
<3 • • • 


a i,i+2 





52) 



and 



^i,i+2 
M+l 



l,i+2 



-j,i+2 -j,i+2 
a i+2,i+l a i+2,i+2 



-j,i+2 -j,i+2 

n n 

"1,2+1 "1,1+2 

^i,i+2 ^',1+2 



"2.i+l 



"2,i+2 



•1)* 



i,*+2 

2 1,2 



Z i,2 



i,*+2 
2 l,i+l 



a 



3,i+2 



2 1,3 



i,*+2 

n,3 



i,i+2 
Z l,i+2 



j,*+2 
l i,i+2 



(83) 



We shall now discuss four specific identities that are going to be very useful in next section. 
The first two identities are relevant to verify the exchange property of the two-particle state under 
the respective rapidities. Both of them are derived by selecting the same indices % = j = 1 in 
Eqs.f l8~Tf82l . Considering that R(X, /J,)i^R(fi, X)i\ = 1 it is possible to write Eq. flHTj) with i — j — 1 

as, 

R(\, fifcl R(ji, X)\j 



while Eq.(J82]) with i 



Dr\x,»)D^\n,X) 
j = 1 gives us, 



R(X, iif 2 '] R{fi, A)a'i 



2,1 ' 



4) 



■^§D^(X^). 



5) 



The two remaining examples are going to be used in the three-particle problem. As before 
they are important to demonstrate the rapidities symmetrization properties of this state. They 
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are derived from Eqs. (182(1831) by choosing the indices i = 2 and j ; = 1, namely 







R(Xt, 




i?(A 1: 


A)|a 






-R(A, Ai)^ R(X, Ai)g| 




A) 


R(Xu 


A)S 


i?(Ai, 


A)^3 


— 




i?(A, Ai)^} i?(A, Ai) 3 | 2 




A) 


R(Xx 


A)S 




A)g 






-R(A, Ai)^ i2(A, Ai) 3 '2 






R(Xi 




i?(Ax. 


A)^2 






-R(A, Ai)^} i?(A, Ai) 3 ' 2 






R(X 1; 


A)g 


i?(Ax. 


A)JJ 




R(X, Ai)^ -R(A, Ai)^ 


d? 2) (Ai 


A) 


R(Xi, 


A)S 


i?(Ai, 


A)S 




-R(A, Ai);^ i?(A, Ai)g'2 


Dr\X h X) 


R(Xi 


A)J3 


i?(Ai, 


A)B 




^(A, Ai)^ i?(A, Ax)^2 






R(Xi 


A)g 


R(Xx, 


A)^2 




-R(A, Ai)^} i?(A, Ai)^2 



i6) 



57) 



In order to construct further relevant identities we found the necessity of combining the results 
of two systems of linear equations. The first of them is obtained by means of the linear combination 
U[b i k + 2]^^-U[b i i + 2]^^ +3 



for b = 1, 2, 3 and k = 1 i — From these equations we 

e,k+2 Uj l,i+2~ Uj e,i+2 u 'l,k+2 lwl ° ~~ 1,1+1,1+2, 



are able to write a 3x3 system of equations for variables a^'l^Oi'*!^— ai'^Oi'tt^ for e 



namely 
/ 



,j,i+2 J,i+2 



7\i+2 i\i+2 



l 2,i 



l 2,i+l 



u l,i+2 

j,«+2 
a 2,i+2 



j,i+2 j,i+2 i,«+2 



"3,i 



"3,1+1 



l 3,i+2 



/ - J,i4 


-2 

h2 a l,i- 


\-2 
f2 


-a- 


-J,»+2 
a i+l,A;- 


|-2 a l,i- 


h2 
f2 


— a 


-i,i+2 
\ a i+2,k- 


^2 a l,i- 


h2 
t-2 


- a 



-j,i+2-j,i+2 
i,i+2 a l,k+2 



i+l,i 



2 a l,fc+2 



i,i+2 -i,i+2 
i+2,i+2 u l,fc+2 



i-1 

E 

e=2 



/ i,*+2 

a l,i+l-e 

j,*+2 
a 2,i+l-e 

i,*+2 
a 3,i+l-e / 



X 



a 



-j,i+2 
a l,i+2 

t+l-e,i+2 



^>-+2 
l,fc+2 



fl i+l- 



-e,fc+2 



+ 4,1 






-j,i+2 
\ °l,i+2 / 



The next system of equations is generated by solving Eq. 
af i+2^i'fe+2- us i n § Cramer's rule we find, 



for the first variable al'^t^u+l ' 



i-l 


.7,1+2 
a l,j+l-e 


i,*+2 
a l,i+l 


J,*+2 
a l,i+2 


E 


i,«+2 
Q 2,i+l-e 


J>i+2 
a 2,i+l 


i,*+2 

a 2,i+2 


e=l 


u 3,i+l-e 


j,i+2 
a 3,i+l 


i,i+2 
a 3,i+2 



^1,1+2 



-i,i+2 
Z l,fc+2 



-j,i+2 -j,i+2 
a i+l-e,i+2 °i+l-e,fc+2 



°k,l a l,i+2 



for A; 



j,i+2 j,i+2 

a l,i+l a l,i+2 

3,i+2 j,i+2 

a 2,i+l a 2,i+2 



i-l, 
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Note that Eqs. (l89l) can be interpreted as a system of linear equations where the variables are 
the 3x3 determinants. In order to obtain our final identity we first solve Eqs. (l89|) for the 3x3 
determinant indexed by e = 1. Next we make the ratio of this solution with Eq. flHTl) considering 
the replacements i — > i + 1 and A <-> \i. After some cumbersome simplifications we find the 
following identity, 







i,i+2 
a l,i+l 


J.i+2 
a l,i+2 






<*2,i 


J,i+2 
°2,i+l 


J.i+2 
a 2,i+2 


a l,i+l 




a ii 
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a 3,i+l 


j.i+2 
a 3,i+2 




7,i+l 


i,i+l 
°l,i+l 




j.i+2 
°l,i+l 


i,*+2 

a l,i+2 


«2,i 


a 2,i+l 




i,*+2 

a 2,i+l 


i,*+2 
a 2,i+2 



<2 



%2 



2 l,i+l 



a- 



Hi.*+2 

l,i+2 
2,i+2 



7,3, i+2 
l,i+l 

^J,i+2 
2,i+l 



a 



a 



-i.i+2 



(7/ 



1,4 
^j,i+2 



2.1 



2 1,3 



-j,i+l 
X l,i+1 



i— 1,3 



i,i+2 
i— l,iH 


-2 


a' 


— 7,i- 


-2 
h2 


a 


«2,i- 


Is3 


a 



i-l,i+l 



-j,i+2 
2 i-l,4 



l,i+l 



2,i+l 



^i,i+2 



1,3 



;=j>*+2 -j,i+2 



%i+2 



%i+l 



2,3 



-j,i+2 

< 3 
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We close this section by presenting the explicit expressions of two other identities in terms 
of the i?-matrix elements. These relations are going to work together with those coming from 
the Yang-Baxter equation in order to simplify the eigenvalue problem. In what follows we shall 
present them in a notation that will be useful for next subsection. The first one follows from 
Eq. (l9"0~]) choosing j = 1, 



(i+1,2). 



)(i+2,4) 



(A,/i) 



,(i+l,3) 



(i+2,3). 



(91) 



where the determinant d\ (A, /j) is given by, 



#0,A)H 



for i < N 



(92) 



■ • • X )b,i+i-b 

The second one follows from Eq. fl82|) also with j = 1 and by performing the rapidity exchange 
A<^/iwe find, 



D 
D 



2 

(i+1,0) 



(A,/i) 

(A,/i) 



(i+2,2) 

5 

(i+2,3) 



(A,A*) 
(A,/i) 



(93) 
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such that the determinant Dg \\, //) is, 



(i+2,2) 



(A,//) 



i?(/.,A)J: +1 i2(,x,A) 



4,i-l 
3,i 



A)™ 



for i + 2 < N. (94) 



4.2 Yang-Baxter equation 



In practice to extract informations from the Yang-Baxter equation one has to project out the N 3 x 
N 3 matrix (j4j) on the Weyl basis. Considering the projection on row [(ai — l)iV 2 + (a,2 — 1)N + 03] 
and column [(ci — l)iV 2 + (c 2 — l)iV + C3] we find, 



J2(A X) A a )^i2(Ax, AaJSj i2(A 2j ^22 = ^ 2 ' ^feS^Ai, A 3 )^ J R(A 1 , A 2 ; 



Cl,C 2 

61,62' 



(95) 



where sum on the repeated indices hi is assumed. 

In order to describe our approach we shall denote Eq. (1951) for a given set of indices a, and q by 
the symbol ^-^'^'^(Ai, A 2 , A3). In general, we are required to elaborated on Eq. fl95|) to obtain 
suitable weights identities that are directly useful in the eigenvalue problem. In fact, there exists 
only two exceptions that follow from Eq. (1951) without further manipulations. These identities come 
from the entries Y B^\(\ x , A 2 , A) and YB N Z^\{\ ^2)- After appropriate normalizations they 
are given by, 



i?(A2, A) 2 .'} R(X\, A 2/ 



3,1 



R(^2, A)^ 2 R(\i, A) 3 ^ -R(Ai, A 2 ) 3 ^ R(Xi, A) 21 
i?(A 2 , A)^i -R(Ai, A)a*i i?(Ai, A 2 ) 3 |l -R(Ai, A)|} 



and 



R{X, A 2 )^2 -R(Ai, A 2 )s'i R(Xi, X 2 )l'^ R(X, A iM ■ :i 



2V,3 



-R(^) ^l)jV-l,3 ^(^) ^s)^,! 1 ' 2 



R(\ A 2 )^'i -R(Ai, A 2 )|i 



i?(A 1 ,A 2 )^i?(A,A 1 ^' 2 



N,2 



R(X, Ai)^o 



-R(A, A 2 )^'i 



(96) 



(97) 



The main feature of Eqs. (l96ll97p is that they factorized two different types of weights products 
into a single product. Identities of this sort are essential to carry out the necessary simplifications 
on the transfer matrix eigenvalue problem. This task, however, will require more complex classes 
of products factorization which involve at least three distinct terms. In order to derive such 
complicated relations we have to make combinations of a special set of equations which follows 
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Equation 


Row index 


Column index 




< fc < m{a- l,N-3} 


< j < a - 1 


YB zVkfi,-k,\ ( A l> A ' ^2) 


< fc < m{a- l,iV-3} 


J = -2 
< j <m{a-2,iV-4} 



Table 7: Two families of linear combination used to derive factorization identities of three weights products. 

from Eq. (|95p . In table ([7]) we have summarized two families of linear combinations that will lead 
to three terms factorization identities. 

Let us now explain how to use the linear combination described in the first row of table (J7|). 
We start by first considering the situation when row index is k = 0. In this case we are left with 
equations YB 2+ Jf~ h2 {Xx, A, A2) (k = 0) for < j < a — 1 which can be combined as follows, 

f R(X 1 ,X)^ hl ... R(Xi, A)|° \ ( R(\,\2)%R(\i,\ 2 )Z 1 ' 2 \ ( <i° ^ 
^ R(Xi, X)lXi,i ■ ■ 



^(A^A)^ 1 ' 1 J \ i?(A,A 2 )";l J R(Ai,A 2 )3^ J \ 1 J 

for 2 < a < N - 2, (9£ 



\2,2 



0,a-l 



The right-hand side components of the nonhomogenous system (19"51) are particular cases of 
more general elements v 2 'i, 

a-j+l 



k,3 



(= J2 R(X U A)3+^ fe2 J R(Ai, A 2 )^ 3 '!^t 2_b2 ^(A, A 2 )b~^ j+2 _ 62 , for * = 0, . . . , a - 1 (99) 
62=1 

The main purpose of the linear combination (19"51) is that it provides us the means to get rid of 
the right-hand side of equation Y B 2+ Jf'~ : ' ,2 (Xi, A, A 2 ). For the remaining values k = 1, . . . , a — 1 
we are therefore asked to eliminate the following terms, 

a— k 



ED/\ \ \b 2 ,a-k+l-b 2 d/ \ \ \a+2-b 2 ,2 \\2+j,a 
rC{A, A 2 J _ A . )1 -K(Ai, A 2 J 3+fc]0 _ fe+1 _ 62 -K( i Ai, AJ a+2 _ 



62,^2 ' 



(100) 
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In order to do that we first rewrite the sum fllOOp in an equivalent matrix form, namely 



a—k 



Ep/\ \ \f>2,a-A:+l-&2 p/ \ \ \a+2-i 



b 2 ,2 
k+l-b 2 



6 2 =1 



2../ \ 



for 1 < k < a - I. (101) 



We next note that the right-hand side column of Eq. 



is exactly the & 2 column of the 



matrix defining the linear system f[98|) . If now we substitute the last column of the a x a matrix 
of Eq. fl98|) by the sum (1 10 II) we will end up with a new matrix whose determinant is null since we 
are computing the sum of a — k determinants with two equal columns. Therefore, to eliminate 
the terms fllOOp for k = 1, . . . ,a — lwe just have to replace the last column of the matrix defining 
the linear system (1981) by the left-hand side of equations Y '-BS^'o-i' i(Ai, A, A2). This leads us to 
the following condition, 



-R(Ai, A)^f x j 



R(^i^)l', a a-i vJJ 

^(Ai,A)3^:i 



In order to have a single relation for all values of k we proceed as follows. For k = we solve 
the linear system (1981) for the last component R(X, A 2 )„ i-R(Ai, A 2 )3i while for 1 < k < a — 1 we 
single out the element t^f from Eq. (11021) . After considering these steps we find, 

a+1 

EDA \\a+3-& 2 >&2 p/\ \ \2,a+2-fc 2 p/\ \ \a,2 _ p qO , cfe 

-K(Ai, AJ 3+fe a _ fc K{M, A 2 ) a+3 „ b2il -K(A, A 2 J b2ia+2 _ 62 - Ofc,o ^2,0 + ^2,1' 

6 2 =1 

(103) 



fc,a— 1 



for 1 < < a — 1 , 



(102) 



for k = 0, . . . , a - 1 and 2 < a < N 



2. 



The dependence of functions 5% and Sfi 011 the -R-matrix amplitudes are 



CO 
°2,0 



-R(A, A 2 )^}-R(Ai, A 2 )g^ 



(a+1,3) 



(A, A, 



(104) 
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and 



°2,1 



#(^1> A)o+l!l ••• R&litytl-l V 2,l 



R(\l,\)l a a -t 



a+1,1 k,a— 1 



£>| a+1 ' 3) (A, Ai 



(105) 



such that the determinant D±' (A, Ai) is defined by Eq.(|92l). 

We now turn our attention to the linear combination associated to the second row of ta- 
ble ([7]). As before we have to make a distinction between the cases k = and that involving 
the remaining values k = 1, . . . , a — 1. In the former situation we have to combine equations 
YB^\ ~ {X\, A, A2) (k = 0) with j = —2 and < j < m{a — 2, N — 4} which can be arranged 



as, 



/ 



2,o+l \ 



#( A l> A )o+2,l ••• R(^l)X)t% 



4,o-l 



( 



R{\ ^2) a 'l-R(Al, A 2 )3^ 



0+2,1 \ 



a+1,1 



^(Ai,A)™y\ i^A^^A^ J j 

for 2<a<iV-2, (106) 



/ ,,0,-' \ 



'1,1 

0,0 



0,a-2 



where the components are obtained from the following more general elements, 



a-j-l 



# = £ ^i^St**^^ for Jfc = 0,... 



a — 1. 



(107) 



62=1 



By the same token, the corresponding algebraic manipulations for the values k = 1, . . . , a — 1 
consist in the substitution of the last column of the matrix (11061) by left-hand side of equations 



YB^kli-k 1 C^i> A> ^2) with j = —2 and < j < a — 2. This leads us to the condition, 



#(Ai,A) 

0+2,1 ••• R(^i^)a'1 



2,a+l fc,-2 



-1 



-R(Ai, A) f 2 ,1 



#(Ai,A) 4 ' a vu 



fe,0 



73/ \ \\a+2,l D/ \ a 

-K(Ai, AJ a+2 1 ... it(Ai, AJ 4 x V lx 



a+2,1 fc,a-2 



for 1 < fc < a — 1 , 



(108) 
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Now by solving Eq. (11061) for the component R(X, A 2 )^' \R(Xi, \2)l'\ and by calculating t>ff from 
Eq. ffT08|) we find, 
0+1 

EDH \\a+3-b 2 ,b2 tj( \ \ \2,a+2-6 2 \ \ \a+l,l _ r cO , ofc 

for fc = 0,...,a-l and 2 < a < N - 2, (109) 



6 2 =i 



where functions S® and are given by 



cO 

°i,o 



;-l) a+1 i2(A l A 2 )Si2(Ai,A 



2J3,1 



^^(A.Aa) 



and 



i?(A 1; A) 

a+2,1 • • • 

i?(Ai,A) 

4,o-l 



-R(Ai, A) 

a+2,1 • • • 4,0-1 



4,0-1 



fe.O 



r>/ \ \\a+2,l d/\ \\H--i,i K,o 



a+2,1 fc,a-2 



(a+2,4) 



(A, Ai 



such that the determinant D^ +2 ' 2 \X, Ai) is defined by Eq. (!94|) . 



(110) 



(HI) 



At this point we have at hand the basic ingredients that are necessary to obtain useful iden- 
tities for the eigenvalue problem. We just have to construct the total number of 'a' 2 x 2 linear 
systems of equations by considering the expressions (I103|ll09p . These a systems are composed by 
Eqs.f TT03lfT09|) for < k < a - 1, namely 



R(X, A 2 )a+i ) i R(X, A 2 )^2 



R(Xi, A 2 )2'i-R(Ai, A) 



\2,a+l 



o-l 



R{Xi, X 2 ) 3 \ 1 R{X lj X) 3+k a _ k 

I + 0fe,0 

6 2 / V ^1,1 / V ^1,0 

for < fc < a - 1 and 2 < a < JV - 2. (112) 



-R(A, A 2 )"+}|} i?(A,A 2 )^ 



a+1,1 



62=1 



cO 
°2,0 



To obtain the desirable identity we proceed as follows. We first compute the products 
R(Xi, A 2 )3'i-R(Ai, A)g^ fc a _ k by applying Cramer's rule in the system (11121) . This will lead us to the 
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following relations 



a 

£ 

6 2 =i 





,a,2 
'a+1,1 


R(\ A 2 )^ 2 a+2 - 


-62 


R(\,\ 2 


4 a+1,1 
'a+1,1 


R(\ A 2 )&* a+2 - 


-62 






R(\ -^2)0+1,1 








R(K ^2)0+^1 


3 



\ \2,a+2-fc 2 r»/ \ 4\a+3-6 

-KlAi, A 2 J a+3 _ 62>1 it(Ai, AJ 3+A . a _ 



+ 4,0 

for < k < a - 1 and 2 < a < N 



R(\ A 2 )„+l'i S% 



2. 



(113) 



The next step is to solve the linear system (11131) for the variable b 2 = a which allows us to 
eliminate function D^'°\\, \2)R(\i, ^s'i- By normalizing the result by D± 2 ' (A, Ai)i?(A, A 2 )^'} 
R(X, A 2 )a+ 1 '}-R(Ai, A 2 )g j and after few manipulations we find that, 

-R(Ai, A 2 )^ 2 



a-1,2 



^ 0) (A,A 2 ) 



g^A^g^^AA) Z^, a+1 ' 2) (A, Ai) -R(A, A 2 )°^ 

i?(A 1 ,A 2 )Jl i?(A,A 2 )^ J D( a+2 ' 3 )(A,A 1 ) L>1° +1 ' 3) (A, Ai) ^A^J 

^(Ai,A2) 2 ; 2 Dr i ' 2 )(A ) A 1 )M aW) (A ! A 1 



+ 



^(Ai,A2)S J Dr i ' 3 )(A,A 1 )^ a+2 ' 3) (A,A 1 



for 2 < a < N — 2. 

(114) 



We finally discuss the special situation when a = N — 1. It turns out that the steps in this 
case are fairly parallel to those already described for 2 < a < N — 2. The respective identity is 
obtained from Eq. (11141) by choosing a = N — 1 except for the terms involving the determinants 



D 



(iv+1,6). 



(A, Ai) and D^ N+1 ' 2 \\, Ai). These determinants are then replaced by the following analytic 



continuations, 



D 



(iv+1,6) 



Dt +2 ' b) (\,\i) 



(A, Ai) = lim 

i2(Ai,A) 



a+2,1 
a+2,1 



-1) 



N+l-b 



-R(Ai, A) ^2 

b,N+2 



N+2-b 
',2 

N+2-b 



-R(Ai, A)^' 2 



#(Ai,A) 



(1 



5) 



42 



and 



D 



(TV+1,2) 



(A, A a 



D (a+2,2) 

lim — 

o->N-l i?(Ai, A) 



(A, Aj 



0+2,1 
a+2,1 



-R(Ai, X)1\n-i 



i?(Ai,A) 



fl(Ai,A) 



4,7V-2 
AT,2 

4.JV-2 
3.N-1 



i?(Ai,A) 



H(Ai,A) 



JV,2 
N,2 

TV, 2 
3JV-1 



(116) 



Though the element i?(Ai,A)^t^ in Eqs.f lll5)116j) are not defined we observe that in the 



identity (11141) such determinants always appear as ratios 



D 



(JV+1,4) 



(A,Ai) 



D 



(N+1,2) 



(A,Ai) 



. This means 



that the common element R(Xi, A) jv+i'i wm ^ e cance l e d out in such ratios. We end by mentioning 
that the very particular case a = 2eiV = 3is obtained by setting D^' 4 \x, Ai) = 1. 



5 The eigenvalue problem 

In this section we shall consider the diagonalization of the transfer matrix (1111) related to the 
inhomogeneous U(l) vertex model defined in section [2j The £7(1) invariance of the respective 
statistical weights implies that the transfer matrix ffTTj) commutes with the total spin operator, 

L 

[T(A),5>*] = 0. (117) 

8=1 

A direct consequence of Eq. (lll7p is that the Hilbert space can be separated in disjoint sectors 
labeled by the eigenvalues of the spin operator, namely 

L 

Y,Sn<5>n) = [Ls-n]\<S> n ), (118) 
i=i 

where |$ n ) denotes the eigenvector on the n-th sector. 

In terms of the monodromy matrix elements the transfer matrix eigenvalue problem can then 
be defined by, 

N 

^r aja (A)|$„) = A„(A)|$ n ). (119) 



0=1 
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Here we are going to solve the eigenvalue problem fl 11 9 [) by solely relying on the three different 
families of commutation rules of section [3], the identities between the i?-matrix elements of section H] 
as well as on the pseudovacuum state properties (112p . To this end we will seek the eigenvectors |$ n ) 
in the form of multiparticle states which can formally be represented by the following expression, 



where the rapidities Ai, . . . , X n parameterize the momenta of the particles. 

For consistency the vector with null number of particles <po = 1 is directly identified to the 
reference state |0). The property ( fl6i) concerning the spin of the monodromy matrix elements 
suggests us that the mathematical structure of the vectors 0„(Ai, . . . , A n ) for general n can be 
searched with the help of linear combination of the product of the creation operators. In what 
follows we are going to argue that such vectors are indeed given in terms of a recurrence relation 
involving the N — 1 creation fields 7^ (A). 

5.1 The one-particle problem 

The one-particle state n = 1 corresponds to an excitation of spin s — 1 over the reference state 
|0). From property (fT6l) we see that among the basis vectors T^A) the operator responsible for 
this excitation is 7l )2 (A). The one-particle vector 0i(Ai) is then given by, 



The solution of the eigenvalue problem f)119p for the one-particle state is achieved by using 
the commutation rules between the operators T a ^ a (X) and 7^2 (A). Taking into account the form 



$n> =0n(Al,...,A n ) |0), 



(120) 



0i(Ai)=T lj2 (A 1 ). 



(121) 
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of Eqs. (l27H29l) as well as the pseudovacuum properties (fT2|) we find, 
Ti,i(A) |$x) 



i2(A 1; A)}^ , -R(Ai, A)£ 2 

W H A )~ — ^2T ~ Wl ( Xl )-^ — ^TT T i>2(A) |0) 



R(Xi, A) 2 \ 



(122) 



l r ll (\)Di a '°\x,X 1 ) + Wl (A0^^%gT a , a+1 (A) |0) 



-R(A, Ai)^ +1 ; x 

"■^■AO ^^tai 7~ a _ 1|a (A) |0) , for 2<a<JV-l, 



-R(A, Ai)"j 



(123) 



T N , N (X) 



A)"TT7rT~TT I^O-WaCAi) ^;^! ?at-i,at(A) |0> . (12 1) 



-^(A ) Ai)j V ' 1 i?(A, Aijjy^ 

It is possible to gather Eqs. (ll22H124p into a single expression for all values of the diagonal 
index a. This requires first the definition of a discrete function that is able to project out a set of 
indices. Keeping in mind the extension to multi-particle states we shall define this function by, 



fc=i 



hJk- 



(125) 



By using identity (1791) to reorder the rapidities of the last term in Eq. (11221) we find that 
Eqs. (I122til24[) can be recasted as, 

%, a (X) |$i) = Wa (A)P a (A,A 1 )|$ 1 )-5fr a , a+1 (A) U ; 1 (A 1 ) 1 ^ 1 (a) (A,A 1 )|0) 

- ^T a _ 1 , a (A)«; 2 (A 1 )o^ 1) (A,A 1 )|0), for 1 < a < N, (126) 

The polynomials P a (X, fi) are the terms proportional to the eigenstate in Eqs . f 1 1 2 2Tj 1 2 4 [) . 
namely 

( Riu.X^A 

for a = 1 



P tt (A,/i) 



R{p, A)5;}' 

^ a,0) (A,/i), for 2 < a < N - 1 

TV, 2 
TV, 2 



(127) 



R(\v)7' 2 



for a = N. 



The auxiliary functions o.T'i'^A, //) an d i-^i (A, /i) are proportional to the terms not parallel 
to the one-particle state. This notation for these off-shell amplitudes have been defined for later 
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convenience. They are given by, 



o^(A^) = -^-W) = fMSS ( 128 ) 



The solution of the eigenvalue problem is obtained by summing Eq. (11261) over the index a. 
This gives us the action of the transfer matrix on the one-particle state, 

N N-l 



(V 

a=l a=l 



T(A) |$!> = ^ Wa (A)P a (A, Ax) + K(A!) - w 2 (Ai)] ^o^r(A, A 1 )T , a+1 (A) |0) . (129) 



We see that the terms 7^ a+1 (A) |0) for a — 1, • • - , N are undesirable states because they are 
not proportional to the one-particle eigenstate |3>i). Note that these states are built up from the 
many possible creation fields possessing azimuthal spin s = 1. All of them are eliminated by 
imposing that the rapidity Ai satisfies the one-particle Bethe equation, 

^ttt = L < 130 > 

w 2 {Xi) 

which implies that the corresponding one-particle eigenvalue is, 

TV 

A 1 (A) = ^)«; a (A)P a (A,A 1 ), (131) 

o=l 

We would like to close this section by commenting on the meaning of the notation we have 
introduced for the off-shell function ^' (A,/i). The two right-hand side indices a and b — a are 
directly related to the companion unwanted operator 7^&(A), making explicit its azimuthal spin 
value s a) b = b — a. The remaining index c accounts for the number of weight u>i(Aj) that is present 

c 

in the respective undesirable term of form TT Wi(Xi)T a)b (X) |0). We shall see that this is indeed a 

8=1 

suitable notation to accommodate similar results for general multiparticle states. 
5.2 The two-particle problem 

The two-particle state lies on the n = 2 sector and is constructed by considering the possible 
excitations with spin s— 2. From Eq. ffToT) we conclude that all products of form / 7i fc 1 (A mi ) / 7^ fc2 (A m2 ) 
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whose indices is constrained by k\ + &2 = 4 can in principle contribute to such state. Therefore, 
the most general ansatz for vector 2 (Ai, A2) one could start with is, 

2 N 

02(Ai,A 2 )= ^2 Y2 4^fe m2 (Ami, A m2 )Ti ifcl (A mi )Ti jfc2 (A m2 ), (132) 

mi,iri2=l fci,fc 2 = l 
m\^m2 fci+fc2=4 

where cj^'J" 2 '' ( A mi , A m2 ) are the coefficients of an arbitrary linear combination. 

However, the commutation rules of the operators T11 (A) and T\${n) as well as between T\ 2 (A) 
and T\ i{}i) shows us that not all the terms entering in Eq. (11321) are linearly independent. In fact, by 
using the commutation relations ( I3TH461) we note that the products Tl ) i(Ai) / 7i ! 3(A 2 ), ^i,i(A2)'7i,3(Ai) 
and / 7i ) 2(A2)7i, 2 (Ai) can be eliminated from the linear combination (11321) . Taking this observation 
into account we conclude that an educated ansatz for the two-particle vector is, 

2 (Ai, A 2 ) = r 1)2 (A 1 )T li2 (A 2 ) + gi 1 \\ 1 , A 2 )T 1 , 3 (A 1 )T 1)1 (A 2 ) + gf\\ x , A 2 )T 1 , 3 (A 2 )T 1 , 1 (A 1 ), (133) 

where g${\\, A 2 ) are arbitrary functions to be determined. 

The next task is to investigate under which conditions the ansatz fll33j) is able to provide us 
the two-particle transfer matrix eigenstates. In other to do that we have to commute the matrix 
elements T a>a (X) with the creation operators present in the two-particle ansatz </> 2 (Ai, A 2 ). Let us 
start by describing the steps necessary to disentangle the action of T a>a (X) on the first term of the 
two-particle ansatz (I133p . In this case we have to use all the three classes of commutation rules 
presented in section 3. The first step is to move the diagonal operators through the creation fields 
Tl 2 (Ai) and Tl^A^ with the help of relations (I27H29I) . This generates products between creation 
fields 7[, 2 (A) and T a _i a {ji) that are not ordered in the unique desirable form 7^_ 1>a (A)7i >2 (Ai) or 
^a-i,a(A) / 7i i2 (A 2 ) for a = 2, . . . , N. This common rapidity ordering are then established by using 
the commutation rules (I38|41l) among the fields 7i )2 (Ai) and 7^_i )(1 (A). The final step concerns with 
the presence of annihilation fields on the middle of products involving three monodromy matrix 
elements that are needed to be carried out to the further right in order to act on the reference 
state |0). This is sorted out by employing the commutation rules (|52|) between the annihilation 
fields with the creation operator 7i i2 (A 2 ). The action of T aA (X) on the second and third terms of 



47 



the two-particle ansatz H133[) is however much simpler. In this case we should not worry about 
the right-hand side diagonal field 7x,x(Aj) since its action on the reference state |0) is known to 
be given by Eq. (112j) . Therefore, the respective computation involves only the set of commutation 
rules among the operators T a ,a{X) and Ti^{\) for i = 1,2. This task is easily accomplished with 
the help of Eqs.flMKi- 

A direct consequence of the above described commutations is the generation of several types 
of terms that are not proportional to none of the terms present in the two-particle ansatz f l 1 3 3 f) . 
Among many different types of such undesirable terms there exists a very special family that need 
to be single out first. This family of unwanted terms are always produced by the same diagonal 
operator T a a (X) for a = 2, . . . , N and carries a dependence either on the couplings g% (Ax, A 2 ) or on 
the specific amplitude ratio ^'^(Ax, X?) / R%\{\\, A 2 ). Therefore, the elimination of such unwanted 
terms are only possible after an appropriate choice of the coefficients ^^(Ai, A 2 ). From now on we 
shall refer to the class of undesirable terms that are responsible for fixing the linear combination 
as easy unwanted terms. In the two-particle problem there are four different terms of this sort 
and their structures are, 



| -^(A) Ax)",! 1 ' 2 
-R(A, Ax)"'x 



iui(A 2 



^(Ai,A a )|| 



SX), 



a-2,3 



Wl(A 2 )w 3 (Ai 



^(A, Ax)a,x 
R {\ Ax)„j 



\2,2 



R^U^J 3.1 + (D (A A 



r a _i, a (A)r 2 , 3 (Ai) 

T a _ 2 , a (A) 



^i(Ai) ^^%x ^ 2) (A 1 ,A 2 )r Q _ 1 , a (A)r 2 , 3 (A 2 



R(X,X 



2)a,l 



a-2,3 



. Wl (A 1 ) W3 (A 2 ) ^; A f"' a 1 1 # (Ax, A 2 )T a _ 2 , a (A) (134) 

K \X, M) a ,l 

From Eq. (11341) we conclude that all the above easy unwanted terms are canceled out provided 
we choose functions g% (Ax, A 2 ) as 



S (1) (Ai,A 2 ) 



R(Xi, X 2 )l' 2 l 
R{Xi, A 2 ) 3 'x 



^(A 1 ,A 2 )=0. 



(135) 



Now, by comparing Eq. (11281) with Eq. (11351) we see that the coupling g^(Xi, A 2 ) is exactly the 
amplitude (Ai, A 2 ) introduced in the one-particle problem. This means that the two-particle 
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vector can be written as, 



2 (Ai, A 2 ) = T 1i2 (Ai)T 1i2 (A 2 ) + A 2 )T 1)3 (A 1 )r 1 , 1 (A 2 



(136) 



The two-particle vector (11361) has the nice feature of being ordered as far as the rapidity Ai is 
concerned. Another important property is its symmetry under the permutation of the rapidities 
Ai and A 2 . This property is derived by substituting the commutation rule (1461) for the operators 
7i,2 (Ai) and 7i )2 (A 2 ) in Eq. (11361) and by considering the definition of iT[ (Ai, A 2 ) given in Eq. (1128|) . 
As a result we find, 



&(Ai,A 2 ) = 0(Ai,A 2 ) 



^i,2(A 2 )'7i >2 (Ai 



D^^(X^> A2' 
D2 '° (Ai, A 2y 



T 1 , 3 (A 2 )T 1 , 1 (A 1 ) 



(137) 



where the exchange function 6*(A,/x) is, 



e{\^) = Df fi \\^) 



R{X,fi)l' 2 2 R(X,n) 3jl 



2:2 



(138) 



R(\, fx)\jR(X, 

The proof of the permutation symmetry is completed by using two identities coming from the 
unitarity property derived in section 4. In fact, if we substitute the identity fl85l) in Eq. (11371) we 
find that 2 (Ai, A 2 ) satisfies the exchange property, 



2 (Ai,A 2 ) = 0(A 1 ,A 2 )0 2 (A 2 ,A 1 ), 



while the second identity fl8~4"l) assures us the consistency condition, 



(139) 



0(Ai,A 2 )0(A a ,AO = l. 



(140) 



The exchange symmetry is very helpful to simplify the expressions for the action of the diag- 
onal operators 7^ a (A) on the two-particle state |$ 2 ). It makes possible to represent the various 
distinct contributions to the unwanted terms, generated by the commutation between T a ^ a (X) and 
02(Ai,A 2 ), in a compact and illuminating form. The details concerning such simplifications are 
rather cumbersome and for this reason they have been summarized in Appendix A. In particular, 
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we show how the identities derived in section 4 are very useful to write a recursive formula for the 
two-particle off-shell amplitudes. In what follows we shall present the main conclusions that can 
be reached from the results of Appendix A. Considering the discrete function (11251) we find that 
T a)a (X) I $2) can be written as, 

2 

r a , a (A)|$ 2 ) = w a (x)Y[p a (x,Xi)\^2) 



i=i 



i, ? =i K {*ii*j)2,l 

- 8^- 1 ' N T a , a+2 (X)w l (X 1 )w 1 (X 2 ) 2 ^ a) (X, A 1; A 2 ) |0) 

2 

- ^T^+^E^WA^ |0) 

i,i'=l 

- ^' 2 T a _ 2 , a (A) W2 (A 1 )«; 2 (A 2 )o^ 2 (a - 2) (A,A 1 ,A 2 )|0), for 1 < a < N. (141) 
where function #<(Aj, Aj) is defined by, 

{9(Xi, AA for z < j 
(142) 
1, for i > j. 

The off-shell functions c-^i ( A ; Ai, A 2 ) for c = 0, 1, 2 represent the three distinct contributions 
of unwanted operators carrying azimuthal spin s = 2. In appendix A we have shown that they 
are given in terms of the following recursive relations, 

TrWfx \ \ \ _ R (^> •r(«+l)r\ \ A 1 ^( A » A l)a+2,l -r(2)/ x x \ / 14 o\ 

0-^2 lA,Ai,A 2 J - — , a+2 ,i0^ 1 (A, A 2 J + — — -— 2ji^i (Ai,A 2 J (14dj 

-K(A, AiJ a+21 K(A, AiJ Q+21 

1^2^ (A, Ai, A 2 ) = o^ a) (A,A 2 ) 1 ^ 1) (A,A 1 ) ^ 2 '^i;l (144) 

^^(A, Ai, A 2 ) = -p^^A, Ai, A 2 ) - ijr 2 a) (A, Xj, Xj) — J I1I177T — r~~^"^<(Aj, Aj). 

i, ; - =1 ^KAj) A iJl,l -n-l A i) A iJ2,l 

(145) 
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The symmetry of the two-particle vector under the exchange of rapidities is expected to be 
reflected in functions ^^(X, Ai, A2) and 23 r( £ l \X, Ai, A2). In fact, from Eq. (1141j) we observe that 
such functions are proportional to products that do not contain function 9(\i, A2). Therefore, the 
invariance of the two-particle state |<&a) by the permutation Ai <-> A2 implies that 



4 a) (A,A 1 ,A 2 )=^(A 1 ,A 2 ) c ^(A,A 2 ,A 1 ) for c = 0,2. 



(a), 



(146) 



We now pause to comment on the results we have so far obtained for the two-particle problem. 
The alert reader will observe that the role of the discrete functions 8£ and 6^ 2 are merely to 
project out undefined monodromy matrix elements. These are operators having either the form 
T a _ iA (X) with a — i < 1 or T a>a+i (X) with a + i > N such that i — 1,2. This remark suggests 
Eq.f ll4ip can still be written in a more compact manner. Indeed, we find that Eq. (1141j) is equivalent 
to the following expression, 



W, 



X 



2 2 m{a-l,t} 

a (X) f[ P a (X, A,) |$ 2 ) E T a-P,a+t- P W 

t=l i=l p=Af{0,a+t-AT} 

i=l,2 \ -r-(a-p) 



l<Jl<-<J (( _p)<2 

)<-<it<2 



i<i(t- P +i 



fc=i 



02-t({Aj}-_^ 5 Jt ) t_ p .F 6 (A, Ajj, . . . , X jt ) 

\ 
J 



i=l 



R(Xi, X jk ) 2 'i 



x 



n ^( A ^) n nn x ,2., 



T 0<(Aj i , Aj 



\ 



J 



t-p t 

n n MAj,,A ifc )io> 
?c=i j=t-p+i 



(147) 



where * in the sum means that any terms with jk = ji for I G {t — p + 1, . . . , t} and k G 

1 i=l,...,n 



{1, . . . , t — p) are different. Moreover, the notation {Ai}j^ V.'", represents that out of the many 
possible rapidities Ai, • • • , A n those indexed by A^, • • • , Xj are not allowed in the set. 

The latter results contains all the basic ingredients we need to solve the two-particle problem. 
The solution of the transfer matrix eigenvalue problem consists basically in taking the sum of 
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either Eg. (11411) or Eg. (1147[) over the diagonal index a. By performing few rearrangements in this 
sum we find that, 



T(A)i$ 2 >=w a (A)n^(A,Ai)|$2> 



i=l 



N-l 2 

EE 

a=l i,j=l 



-^^ A ^ + WA,)o^ ) (A,A l ) 



R(\j, Ai) 2 '! 



^(A», Aj)^i , 



x ^ < (A J ,A i )T aia+1 (A)0 1 (A i )|O) 



7V-2 

E 

a=l 



(Ai)^i(A 2 ) 2 J- 2 (a) (A, A 1; A 2 ) + w 1 (X i )w 2 (X j ) 1 ^ a \\ K *j)0<&, A 



(a)/ 



+ W 2(AiV2(A 2 )o^ a) (A,A 1 ,A 2 ) T a , a+2 (A)|0) 



'148) 



From Eq. (11481) we can see the possible classes of unwanted terms that should be canceled out. 
The first term come as direct extension of the structure present in the off-shell one-particle state. 
The second term is inherent to the two-particle state and corresponds to the three possible manners 
of constructing unwanted terms with azimuthal spin s = 2 of type Zi,a+2(A)wi(Ai)w;j(A 2 ) |0) for 
i, j = 1,2. In order to cancel all such unwanted terms we have to use the properties of the auxiliary 
functions (A, Ai) and 2 -7 r 2 (A, Ai, A 2 ) given by Eqs. (jl28| 11451) . More specifically, considering 



that !^ o) (A,Ai 



F[ a \\\x) as well as by substituting the amplitude 2 JF 2 (A, Ai,A 2 ) in 



(a), 
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Eq. ffT48|) we find, 

2 

T(A)|$ 2 ) = w a (\)l[P a (\,\ t )\<5> 2 ) 



i=l 



N-l 2 



+ EE 

a=l i,j=i 



-R(Aj, Aj) 2 'i 



-R(Aj, A^a'i 



oJ^i ^(A, Aj 



x e<(A i ,A i )T aia+1 (A)0 1 (A i )|O) 



N-2 



+ {K( A l) W l( A 2) -W 2 (Ai)w2(A2)]o-^2 a) ( A > A l> A 2 



o=l 
2 



+ y^wi(Aj 



i ,j = 1 



i?(Aj, Aj) Aj) j'| 

Wl ( A i) x xV.l - WAj)— — — 2T^(A i; Aj) 



-R(Ai, A,-) 2 'i 



-R(Aj, Aj) 2 'i 



-R(Aj, Aj)^! 
-R(Aj, Ai)}|} 



x ^(A, A,, A i )0 < (A i , A,) f T a , a+2 (A) |0) . 



(149) 



From Eg. (1149[) we see that the unwanted terms either proportional to 7^ a+1 (A)0i(A./) |0) or 
to T a<a+ o(X) |0) vanish provided that the rapidities Ai and A 2 satisfy the following Bethe ansatz 
equation.. 



W 2 (\j) 



n^A, 



-R(Aj, Aj)i^ R(\i, Aj)^i 
-R(Aj, Ai)ai -R(Aj, A^K 



i=l 



for j = 1,2. 



(150) 



and consequently the two-particle eigenvalue is, 

N 2 

A 2 (A) =J]«; a (A)JJP a (A,A i ). (151) 

a=l i=l 

The typical feature of integrable theories is that the two-particle results already contain the 
main flavour about the structure of the spectrum. This means that expressions (11511) and (I150p are 
believed to be valid for general n-particle states. To benefit from the knowledge of the eigenvectors 
one has, however, to perform the complete analysis for few other multiparticle states with n > 2. 



( \ \ 

Note that Eq. (|150p implies the relation 1 f — — ^— = 1 thanks to the property 0(Ai, A2)#(A2, Ai) = 1 

±J : W2(A<) 

2—1 V ' 
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5.3 The three-particle problem 

The three-particle state is constructed in terms of a linear combination of products of creation 
fields with spin s — 3. Thus, the most general ansatz for the three-particle vector is, 

3 N 

03(^i> A 2 , A 3 ) = ^ ^ c^' 2 ™ 2 3 ' m3 ^(A mi , A m2 , \ m3 )Ti kl {\ mi )Ti )k2 (\ m2 )Ti^ z {\ mz ), 

m\ ,m,2 ,7713 = 1 ki,k2,k^=l 

(152) 

where 4™fc 2 ™3 (^"h> A m2 , A ms ) are arbitrary coefficients. 

This general ansatz can be further simplified by considering the following two step pro- 
cedure. First we use the commutation rules between the fields Ti^(X) and 7i,b(/i) for b = 
2,..., 4 to eliminate from Eq. (1152ft the linear dependent products 7i i i(A m j7i ) i(A m2 )7i ) 4(A m3 ), 
^,i(Ami)'7i,4(A m2 )'7i i i(A m 3), 7i ) i(A mi )7i ) 2(A m2 )7i ) 3(A m3 ), 7i,i(A mi )7i j 3(A m2 )7i ) 2(A m3 ), Ti^{X mi ) 
x ^i,i(A m2 ) / 7^ i 3(A m3 ) and 7i !3 (A TO j7i ) i(A m2 )7i )2 (A n j, 3 ). The next step is to reorder products of cre- 
ation operators associated to the combinations 7i 2 (A) with 7^2 (jti), T X2 (A) with 7j 3 (/x) and 7^ (A) 
with T12 (/x). This task is performed with the help of the commutation rules given by Eq.(T4T]). 
Taking into account these steps we find that the ansatz for the three-particle vector fl 1 5 2 [) becomes, 

3 (Ai, A 2 , A 3 ) = T li2 (A 1 )0 2 (A 2 , A 3 ) + ^(Ai, A 2 , A 3 )T li 3(A 1 )0 1 (A 2 )T 1 , 1 (A 3 ) 
+ gf\K A 2 , A 3 )T li3 (A 1 )0 1 (A 3 )T lil (A 2 ) + gf\\ Xt A 2 , A 3 )T li3 (A 2 )0 1 (A 3 )T lil (A 1 ) 
+ #(Ai, A 2 , A 3 )r li4 (A 1 )r lil (A 2 )r lil (A 3 ) + gf\\ u A 2 , A 3 )T m (A 2 )T 1 , 1 (A 1 )T 1 , 1 (A 3 ) 
+ 4 6) (A!, A 2 , A 3 )r M (A 3 )r 1|1 (Ai)r 1|1 (A 2 ) + ^ 7) (A 1; A 2 , A 3 )T 1 , 2 (A 1 )T 1 ,3(A 2 )T 1i1 (A 3 ) 
+ gf\\ x , A 2 , A 3 )r li2 (A 1 )r li3 (A 3 )r 1)1 (A 2 ) + ^ 9) (A 1; A 2 , A 3 )T 1 , 2 (A 2 )r li3 (A 3 )r 1 , 1 (A 1 ), 

(153) 

where g 3 ^(Ai, A 2 , A 3 ) are going to be fixed bellow. 

Note that the three-particle ansatz has been written in terms of the one-particle and two- 
particle vectors. In this manner it will be able to use the previous results already established in 
sections 15 . 1 1 and I5T21 Further progress is made by investigating the result of the commutation of the 
diagonal operators with the four distinct types of product of operators defining the three-particle 
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vector (11531) . As in the two-particle problem there is no need to worry about the right-hand side 
diagonal fields 7^1 (Aj). In what follows we describe the main steps one has to perform in order to 
solve this problem. 

Let us start by considering the commutation of 7^ a (A) with the most involved term 7i i2 (Ai) 
x</> 2 (A 2 , A 3 ). The commutation of the diagonal fields with 7^ 2 (Ai) is once again done by using 
Eqs. ( 12711291) . This step generates besides the diagonal fields several annihilation operators acting 
on the two-particle vector 02(A2,A 3 ). More specifically, we have carried on operators of type 
T d+a>a {\) or T d+a a (\i) for d = 0, • • • ,N — a through the two-particle vector 02(A2, A 3 ). This task 
is immediate for d > 3 but it requires a number of computations for d = 0, 1,2. In particular, 
we have to reorder, as far as the rapidity A is concerned, the following products of operators 
7"i,2(Ai)T aiiai+ i(A)Ti i 2(A i ) with j = 2,3 and Ti i2 (Ai)T a2ja2+2 (A) as well as 7" lj3 (Ai)7; iiai+1 (A) where 
dj = a, . . . , a — %. These computations are rather cumbersome and have been summarized in 
Appendix B. 

We now discuss the commutation of the diagonal fields with the second term 7i ) 3(A l -)0i(Aj). In 
this case we start by commuting the operators 7^ a (A) with 7i )3 (Ai) by means of Eqs. (13111341) and 
as before we produce diagonal and annihilation fields. These operators are then carried out to 
the right with the help of Eqs. (l52|ll26p . The final step consists of reordering products of creation 
fields such as 7i j3 (Ai)7^ a+ i(A) and 7i j3 (Ai)7^_i )a (A) having the rapidity A in the further right. 
This problem is resolved by using the commutation rules given by Eq. (JJTJ) as well as by the linear 
system of equations (I42f43l) with d\ = 0, b\ — 3. The commutations between T aa {X) and the 
third operator 7^4 (Aj) are obtained more directly. In the cases a = 1 and a = N they follow from 
Eqs. ( I19|l20l) with 6 = 4, respectively. For 2 < a < N — 1 we have to solve the linear system of 
equations (121l24j) with b = 4 by applying the standard Cramer's rule. 

We finally discuss the commutation of the diagonal fields with the last term 7^ 2 (Aj)7^ 3 (Aj) 
for i,j = 1,2,3 such that i ^ j This turns out to be the simplest case in our analysis due to 
the following reason. The commutation of 7^ a (A) with 7i j2 (Aj)7i i3 (Aj), by employing Eqs. (I27H291 
[3TH34T) . produces a special class of easy unwanted products that are not generated by any of the 
previous three terms entering the three-particle ansatz. The form of such undesirable products is 
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given by, 

• ^ 4+i+i) (A 1; A 2 , A 3 )T a _3,a-i(A)T 2 ,3(A l )T3 i3 (A i ) for a > 3 (154) 

• ^ 4+i+i) (Ai, A 2 , A3)r a -4,a-i(A)T3,3(A,)r 3 ,3(A i ) for a > 4. (155) 

Therefore we must project out the last term Ti2(Xi)Ti^(Xj) from the linear combination (11531) 
which is achieved by setting #3 (Ai, A 2 , A3) = #3 (Ai, A 2 , A3) = gf\Xi, A 2 , A 3 ) = 0. This choice 
prompts a sequence of similar cancellations of other terms in the linear combination (11531) . The 
first one concerns with the product 7i,3(A 2 )0i(A3) which generates the following easy unwanted 
terms, 

• ^ 3) (Ai, A 2 , A 3 )r a _ lia (A)T 2i 3(A 2 )0 1 (A 3 ) for a > 1 (156) 

• #(Ai, A 2 , A 3 )r a _ 2ia (A)r 2i 3(A 2 )r 2i2 (A 3 ) for a > 2 (157) 

• ^ 3) (A 1 ,A 2 ,A 3 )T a _3,a(A)T3 i 3(A 2 )T 2i2 (A 3 ) for a > 3. (158) 

which are eliminated provided that #3 (Ai, A 2 , A3) = 

The next case is associated to the operator Ti^(Xj) for j = 2,3. Considering the above con- 
straints, the commutation of T aA (Xj) with Ti^(Xj) generates the following types of easy unwanted 
terms, 

• <?f +i) (Ai, A 2 , A 3 )r a _ 1)a (A)T 2i4 (A i ) for a > 1 (159) 

• gt^iXi, A 2 , A 3 )T a _ 2 , a (A)T 3i 4(A J ) for a > 2 (160) 
. ^ 3+i) (A 1; A 2 , A 3 )T a _3, a (A)T 4i 4(A J ) for a > 3, (161) 

• ^ 3+i) (A 1; A 2 , A 3 )r a _ 1)a+1 (A)T 2)3 (A j ) for 1< a < N (162) 

• 33 3+j) (Ai, A 2 , A 3 )T a _ 2 , a+1 (A)r 3 ,3(A i ) for 2 < a < N. (163) 

which are canceled out provided we choose gf\X\, A 2 , A3) = gf\X\, A 2 , A3) = 0. 

The restrictions we have found so far bring a considerable simplification to the structure of 
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the three-particle ansatz (11531) . We are now left with only four independent terms, 



03(Ai,A 2 ,A 3 ) = Ti i2 (Ai)02(A 2 , A 3 ) 

+ ^(Ai, A 2 , A 3 )r 1) 3(A 1 )0 1 (A 2 )r 1)1 (A 3 ) + <fe W (Ai, A 2 , A 3 )r 1 ,3(A 1 )0 1 (A 3 )r 1 , 1 (A 2 ) 
+ ^ 4) (A 1 ,A 2 ,A 3 )T m (A 1 )T 1i1 (A 2 )T 1i1 (A 3 ) (164) 



The structure of the three-particle vector (11641) shows some remarkable characteristics. We 
first note that the left-hand side of the linear combination is always commanded by the creation 
fields T\ ,&(Ai) where the index h is limited by the maximum possible spin value. Next, each term of 
the linear combination is multiplied by the previous admissible (4-b)-particle state with rapidity 
A 2 and A 3 . Finally, diagonal fields 7i,i(Ai) are added to the further right to complete the total 
number of three possible rapidities. Altogether, these are rather illuminating features that will 
be of great help to set up the ansatz for arbitrary multiparticle states. To determine completely 
the linear combination (I164p we have to return to the analysis of the easy unwanted terms. Now 
the same form of a given remaining easy unwanted term is generated by two distinct types of 
operators present in the three-particle vector (I164p . In particular, the operators / 71 2 (Ai)0 2 (A 2 , A 3 ) 
and 7l, 3 (Ai)0i(A./) generate the following class of easy unwanted terms, 



a-1,2 



xr a _ lia (A)T 2i3 (A 1 )0 1 (A i )T lil (A i ) for i,j = 2,3;i^ j. 



R(Xi, Xj 



3)2,1 e=2 



(165) 



On the other hand operators / 7i, 2 (Ai)^ 2 (A 2 , A 3 ) and ^^(Ai) are able to produce yet another 
family of easy unwanted terms, 



a-1,2 



3 jf\X u A 2 , A 3 ) - ^(A a , A 2 , A 3 ) T a _ lia (A)T 2i4 (A 1 )T lil (A 2 )r lil (A 2 ) for a > 1. 

(166) 



The cancellation of these two families of easy unwanted terms is obtained by choosing functions 
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9 3 (Ai, A 2 , A 3 ), g ( 3 J (Ai, A 2 , A 3 ), and g ( 3 J (Ai, A 2 , A 3 ) as, 



^(Ai,A 2 ,A 3 ) 



\T \ (Ai, A3) 




fl(A 2 ,A 3 ) 



(167) 



<?nAi,A 2 ,A 3 ) 



*T(Ai,A 2 ) 



-R(A 2 , A 3 ) 2 ^ 
-R(A3, A 2 )i^i 



(168) 




C j (Ai,A 2 ,A 3 ) 



2 (Ai, A 2 , A3). 



(169) 



It is important to emphasize that the choices (I167H169]) eliminate all possible easy unwanted 
terms generated by the commutation of the diagonal field 7^ a (A) with the simplified three-particle 
vector (I164p . By substituting the constraints (1167111691) in Eq. (1164[) we find the three-particle vector 
is now fixed by the expression, 

03(Ai,A 2 ,A 3 ) = T 1;2 (Ai)0 2 (A 2 , A 3 ) 



We now turn our attention to investigate the exchange properties of the three-particle vector 
(11701) concerning the rapidity permutations A 2 <-> A3 and Ai <-> A 2 . In what follows we shall show 
that this vector satisfies the relations, 



The permutation for the rapidities A 2 and A3 follows immediately from the two-particle ex- 
change symmetry. In this case we just have to use Eqs. (ll39|ll46l) . The technical steps to show the 
permutation for the rapidities Ai and A 2 are more intricate. We start by commuting the operators 
7j 2 (Ai) and 7j 2 (A 2 ) in the product '7i j2 (Ai)0 2 (A 2 , A 3 ) by using the two-particle expression (11361) 
and its symmetry property (11401) . We next carry on the operators ^^(Ai) and / T L1 (A 2 ) through 
the operator Ti^{\3) by using Eq. (l27|) . By performing these two steps procedure we find the 



■1,, j— ^ * -j ' 

+ 2-^ 2 (2) (A 1 , a 2 , A 3 )r li4 (A 1 )r lil (A 2 )r lil (A 3 ). 




(170) 



03(Ai, A 2 , A 3 ) = 0(A 2 , A 3 )0 3 (Ai, A 3 , A 2 ) = 9(Xi, A 2 )0 3 (A 2 , A x , A 3 ) 



(171) 
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expression, 



T 1 , 2 (A 1 )0 2 (A 2 , A3) = 0(Ai, A 2 )T 1 , 2 (A 2 )0 2 (A 1 , A 3 ) + 9(X 1 , A 2 )i F? ] (X 2 , A 1 )T 1 , 3 (A 2 ) 
x P 1 (A 1 ,A3)r 1)2 (A 3 )r 1)1 (A 1 )- 1 ^ 1 (1) (A 1 ,A3)r 1)2 (A 1 )r ljl (A3) 

i^i (2) (Ai, A 2 )T li3 (A 1 ) [Px(A 2 , A 3 )T li2 (A 3 )T lil (A 2 ) - 1 J\ (1) (A 2 , A 3 )T lj2 (A 2 )T l5l (A 3 
+ i^i (2) (A 2 , A 3 )T 1 , 2 (A 1 )T li3 (A 2 )T 1 , 1 (A 3 ) 
- 0(Ai, A 2 )i^ r 1 (2) (Ai, A 3 )T 1 , 2 (A 2 )r li3 (A 1 )r ljl (A 3 ). 



(172) 



The next task is concerned with the elimination of products / 7i, 2 (A 2 )7^ 3 (Ai)T L1 (A 3 ) and 7l i2 (Ai) 
x7i i3 (A 2 )7i i(A 3 ) from Eg. ( 11721) . This is accomplished by using the commutation rule between 
7i ;3 (A) and 7i >2 (/i) given by Eq.( l4TT) with d\ = and b± = 3. After commuting these operators in 
Eq.([I72D we find, 

T 1 , 2 (A 1 )0 2 (A 2 , A 3 ) = 0(Ai, A 2 ) fT 1 , 2 (A 2 )0 2 (A 1 , A 3 ) + 1 ^ (2) (A 2 , Ai)Pi(A 1; A 3 )T li3 (A 2 )T li2 (A 3 ) 



x r ltl {\ x ) + X H { ?{X 2 , A 1; A 3 |2)T 1 , 3 (A 2 )T 1 , 2 (A 1 )T lil (A 3 ) + 2 ^>(X 2 , A 1; A 3 )T li4 (A 



(2), 



x T M (A 1 )T 1)1 (A 2 



1 ^ 1 (2) (A 1 ,A 2 )P 1 (A 2 ,A 3 )T li3 (A 1 )r li2 (A 3 )T lil (A 2 ) 



+ i^i 2) (Ai, A 2 , A 3 |2)T 1 , 3 (A 1 )T 1 , 2 (A 2 )T 1 , 1 (A 3 ) + 2 ^(A 1; A 2 , A 3 )T M (A 1 )T 1 , 1 (A 2 )T 1 , 1 (A 



(2), 



(173) 



where functions ihtf\\ Ai, A 2 |2) and ^^(A, Ai, A 2 ) are given by 



(2), 



{H?\\\ u \ 2 \2) 



2^2 {\ ^2) 



^ 2 ' 0) (A, Ai) R{\, A 2 ) 2 ; 2 R(X, Ai) 2 ;} R(X, A^ 2 ; 2 R(X 1 , X 2 )\j 



(174) 



Df ' 0) (A, Ax) P(A, A 2 )J} P(A, Aa)|;J R(X, X 1 f 3 j R(X 1 , X 2 ) 2 2 j 
RiXuX^ljD^iXuX) 

Df 0) (A, Ai) P(A, A 2 ) 2 ; 2 P(A, AQ 2 ;! P(A, AQjg _ R(\ u A 2 ) 2 ; 2 gg^Ai, A) 

•K(Ai, A 2 )^i ^'^(Ax, A) 
(175) 



£>f 0) (A, AO P(A, A a )5J P(A, Ai)£ P(A, A x 



3,1 



\3,1 



'4,1 



We now reached a point in which we have to implement simplifications on the expressions 
(11 7411 1751) . This is done by taking into account identities coming from the unitarity relation. We 
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first use Eg. (1791) to change the order of the rapidities Ai and A2 on the second term of Eq. (I174j) . 
We then consider the results (j86|87p to reorder the rapidities A and Ai on the determinant ratios 

<f Eqs. (117401751 ). After h 
identities for functions {H^ (A, Ai, A 2 |2) and 2 .F 2 (A, ^2) 



D (3,o)^ Al,A ' > an d ^.o)^ 1 '^ °f Eqs . (1 1 74 II 1 751) . After implementing these steps we find the following 



iRf\\ Ai, A 2 |2) = ^i 2) (A, A 1; A 2 |2) = 9(X 1 , X 2 ) ^ U ^^ i^ (2) (A, A 2 ) (176) 
2 ^ 2 (2) (A,Ai,A 2 ) = 2 ^ 2 (2) (A,A 1; A 2 ) (177) 

where iHi(X, Ai, A 2 |2) has been defined in Eqs.l 1A.12l) . 

Now if we substitute the above expressions in Eq. (jl73j) we find that the product 7i >2 (Ai) 
02 (A 2 , A 3 ) can be finally be written as, 

T 1 , 2 (A 1 )0 2 (A 2 , A 3 ) = 0(A 1; A 2 ) [t 1 , 2 (A 2 )0 2 (A 1 , A 3 ) + ijf\x 2 , X 1 )P 1 {\ 1 , A 3 )T 1>3 (A 2 )T li2 (A 3 ) 

xT 1)1 (A 1 ) + 0(A ls A 3 )i^; 2) (A2, A 3 )Pi(A 3 , Ai)r li3 (A 2 )r li2 (A 1 )r 1)1 (A 3 ) + 2 ^ 2) (A 2 , A x , A 3 ) 

x r M (A 2 )T lil (A 1 )T 1 , 1 (A 3 )] - [i^ 1 (2) (A 1 ,A 2 )P 1 (A 2 ,A3)T li 3(A 1 )T li2 (A 3 )T 1 , 1 (A 2 ) 

+9(X 2 , A 3 )i^ 2) (Ai, A 3 )Pi(A 3 , A 2 )T lj g(A 1 )T 1 , 2 (A 2 )T 1 , 1 (A 3 ) + 2 ^ 2 (2) (Ai, A 2 , A 3 ) 

x r M (A 1 )T 1 , 1 (A 2 )T 1 , 1 (A 3 )], (178) 

which is indeed equivalent to exchange property (I17ip . 

We now turn our attention to the action of diagonal fields 7^ (A) on the three-particle state 
I $ 3 ). The permutation property (11711) is once again of fundamental help to gather the contributions 
of T a ^{X) I $3) in closed forms. In Appendix B we have described main technical details of these 
computations and our final results are given by Eq. (1B.16l) . It turns out that this expression can 
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also be rewritten in a neat compact form, namely 

3 



3 m{a— l,t} 

T a , a (A)|$ 3 > = ™ a (A)n^(A,A0|<t>3>-2] ^ T a _ p , a+t _ p (A) 

i=l *=1 p=M{0,a+t-AT} 

^ 03-t({Ai}-^ 2 ' 3 it ) <- P ^ (a p) (A, X jl7 . . . , A 3 - t ) 

l<il<-<J( t _p)<3 

l<J(t- P +i)<-<Jt<3 



X 



X 



X 



t-p 



n 



2T y <l A " A 3k) 



\k=l 
( t 



;=1 ^(Ai, AjJ^ 
i¥=3i,—Jt 



V 



n ^(a,) n 



n IHw ^^in n m^aji°>> 

i=i K V A JP A iJ2,l / fc=n= t _ p+ i 



(179) 



As expected the three-particle state generates new unwanted terms proportional to the creation 
operators T aya+ ^(X) carrying spin s = 3. The corresponding off-shell amplitudes 
C ^-3^(A, Ai, A2, A3) for c = 0, 1, 2, 3 are obtained from the following relations, 



n=3 



>- >- ^ " ^ Al)a ;; .F 2 (a+1) (A, A 2 , A 3 ) + ]T M^ o^CA, A, 



o-T 7 3 (A, Ai, A2, A 3 ) 



n{A, Aij a+31 ,^. =2 -n^A, Aij a+3)1 

X ^ (Ai, A7J— — 2 i y<^Ai, Ajj + +3>1 
^mA, Ajj 2 ,i ^A, AiJ a+3)1 



-R(A«, Aj)H 



-R(A, Ai)^ 3)1 



r2^ 7 2 2 ' ) (^i) A 2, A 3 ) 



(180) 



i2(A 2 ,Ai)i;ii2(A 3 ,A 1 )i;i 



iJF 3 •'(A, Ai, A 2 , A3) — 0^2 ^(A, ^2? A3) + ; (A, Ai 



7 (a+2) 



1,1 



-R(A 2 , Ai)^ -R(A 3 , Ai)^J 



2,1 



(181) 



T {a) ( , XXX) 77(a) a \ \ 77(»+l)M A A x ^(^3 ; Aj)^ -R(A 3 , Ag)^ 

2-^3 IA, Ai,A 2 ,A 3 J - J- 1 [A, A 3 J 2 7- 2 (A, Ai,A 2 j — — 7—^—77 7 TTT 

-n-(A 3 , AiJ 21 rqA 3 > A 2 j2,i 



(182) 



61 



J~s (A> Ai, ^2, A3) — —0^3 ; (A, Al, A2, A3) 

3 2 
- ^ ^ i^3 a) (Aj A J1 ,A; 1 ,A /2 )]^[6' < (A J1 ,A; 



(a), 



^(Aii, ^1)2,1 



2,1 



-R(Aji, A;j|} -R(Az fc , AjJ^i 



1,1 



il = l l<Ii<!2<3 



k=l 



2 , 2 / 2J 3 (A A JD A i 2? A «J ll y <l A i fe) 



1,1 ' 



i<ii<i2<3 «i=i 



k=l 



R( X 3h> X h)2,l -R(Aii) ^j k )l,l 



(183) 



Direct inspection of Eqs. (jl80til83p reveals us that all the three-particle off-shell amplitudes are 
determined by a recurrence relation that depends on the previously calculated off-shell functions 
associated to the one-particle and two-particle problems. From Eg. (11791) we note that the only 
unwanted terms not carrying function 0(Aj, Xj) are those proportional to qJ 7 ^ (X, Xi , A2 , A3) and 
3^3^ (A, Ai, A 2 , A 3 ). Therefore, the three-particle vector permutation property implies the following 
relation, 



c ^ a) (A, Ai, A 2 , A 3 ) = 0(Ai, A 2 ) c ^ aj (A, A 2 , A 1; A 3 ) = 0(A 2 , A 3 ) C ^ 3 W (A, A 1; A 3 , A 2 ), for c = 0, 3. 

(184) 

In order to solve the three-particle problem we have first to perform the sum over the diagonal 
indices a = 1, . . . , iV in Eq. (1179j) . After few rearrangements of the sums entering in Eq. (11791) we 
obtain, 

3 3 N—t t 

T(A) |$ 3 > = w a (X) n P (A, A,) |$ 3 ) - £ E T ^ X ) E E* & M^/fZr ) 



i=l 



t=l a=l 



p=0 l<Ji<-<jp<3 

i<i(p+i)<---<j't<3 



x p-^ r t ^(A) Aj i; . . . , Xj t 
( 

V 



n-i(Aij n " /J i>.( a,. a ; 



X 



=1 -R(Ai, A Js ) 2 ^ 

-R(A>,Aj)H 



/ 



n W2 ( a >) n ^ttw 9 <(v^ 

r=p+l i=i K K A jr» A iJ2,l 



/ 



n n ^<(v>a,j. 

s=l r=p+l 



;i85) 



The next step consists of collecting together the contributions proportional to the same type of 
unwanted terms T a a +e(X). This task is performed by plugging in Eq. (11851) the off-shell amplitudes 
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(A, AjJ, 2^2 L \^, Xj 2 ) and 3J3 (A, Aj 1; Aj 2 , Aj 3 ) denned by Eqs. (ll28fl45lll83p . After some 
simplifications we find the following expression, 

N 3 3 N—t t-1 



a=l 



i=l 



t=l a=l 



T(\) i<& 3 > = E w »w II P *( A > A *) i $ 3) -EE E E 

(p 



i=l,2,3 



p=0 :• .•■ s 3 

l<J(p+l)<-<J*<3 



x p ^(A,A J1 ,...,A,)|n^(A J j n ^¥U e <(^^ 



s=l 



t p 3 

n \\ e <^K) n ^<( a - a 

r=p+l s=l 



=1 -R(Aj, AjJ 2 ^ 

i&u—Jt 

\ 



J 



JrJ 



tyjlt—jt 

t 3 



r=p+l 



i=1 -R(A Jr , Aj) 2 ^ 

i&U—Jt 



n a , x v , tt 1 r ^(Aj, A^)^ -A- R(Xj s , Xj r )i i R(Xj r , AjJ^i 

o(K>K)- 11 w i( A >) 11 PM rw 11 



>i,i 



X 



s=l 



r=p+l 



=1 -R(Aj, A Jr ) 2 | 1 g=1 R(X js , Xj r ) 2 li R{Xj r , XjJ^ 

i&u—dt 



(186) 



From Eq. (jl86p we conclude that all unwanted terms proportional to T aA+t (X) can be eliminated 
by imposing that functions inside the brackets are null, namely 



-pit 



n 



ri n ^A ^^^^ f[e(x jr ,x 



r=p+l i—l 

tyhi—Jt 



R{Xj r , Aj) 2 'i R(Xi, AjJ^ s=1 



3s J 



R(^3r)^3s)l,l R (^3sl^3r)l',l r , 100 A n +1 /1C7\ 

x — ; 2T yt tor t = l,2,3 and p = 0, 1. (187) 



^(^V) ^js^'l Rftjsl <\?"r)ljl 

It turns out that this condition is fulfilled provided that the rapidities Aj satisfy the following 
Bethe ansatz equations, 



Wi(Xj) 



R(Xj, Xi) u R(Xi, Xj) 2 i 



M*j) ly^'^RiX^X^lRiX^X 

i¥=3 



11 
•3 J 11 



for J = 1,2,3. 



We observe that Eq. (jl88j) follows directly from Eq. (11871) for t = 1. To show that the other 
unwanted terms proportional to 7^ a+2 (A) and 7^ a+3 (A) are also canceled out we just have to 
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substitute Eg. (11881) in the left-hand side of Eg. (11871) . As a result we obtain, 

n f[o(x jr A i ) R ^ jr,x ^i R ^ i,Xjr) ^ = n n ^ r ^) R{Xjr,Xi) l R{Xi,Xjr) ^ 

s= i ^(A>, AjJ 2 |i -R(Aj s , Aj,.)^ 

By carrying out few simplifications on Eg. (11891) one is able to show that it is eguivalent to the 
expression, 

t 

J] ^(A > ,A Ji )^(A Ji ,A > ) = l for t = 2,3. (190) 

i,r— p+ 1 

i<r 

which is satisfied thanks to the property (11401) of function 8(X,fi). As a conseguence of the 
cancellation of all unwanted terms we find that the three-particle eigenvalue is given by, 

N 3 

A 3 (A)=^w a (A)nP a (A I A j ). (191) 

a=l i=l 

We conclude by observing that our results for the two-particle and three-particle states are 
already capable to unveil us a common pattern not only for the on-shell properties but also to 
the structure of the off-shell Bethe vectors. It should be emphasized that this fact has become 
possible only because we have been working without using any specific U(l) i?-matrix. 

5.4 The multi-particle state 

The previous solution of the two-particle and three-particle problems provides us basic guidelines 
that are useful to construct multi-particle states. The first lesson we learn is that among the 
possible product of creation operators entering the linear combination there exists a large subset 
that in practice do not contribute to the respective state. These terms produce unigue easy 
unwanted terms that are canceled only by setting the corresponding linear combination coefficient 
to zero. Our results indicate that the surviving product of creation operators possess the following 
characteristics. From Egs. (ll36|ll70p we observe that the left-hand side of such products is always 
governed by the creation fields T^Ai) with 2 < b < m{n+l, N} where n denotes n-particle sector. 
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In addition, each operator Tl^Ai) is then multiplied by the previously determined (n + 1 — b)- 
particle vector such that the spin of the n-particle state is reproduced. Finally, a number of (b — 2) 
diagonal fields Ti t i(\i) are added to the further right to complete the total number of rapidities 
Ai, • ■ • , A n present in the n-particle vector. Combining these features together we find that such 
building blocks can be written as, 

e-1 

Ti^iX^n-tttX^-^J J] for e = 1, . . . , m{n, N - 1}. (192) 

fci=i 

where the indices ji for / = 1, • • • , e — 1 take values on the interval 2 < ji < n subjected to the 
condition 2 < ji < j 2 < ■ • ■ < js-i < n. 

Considering the above discussion n (Ai, . . . , A n ) is obtained by taking the most general linear 
combination of the terms (11921) . namely 

m(n,N-l) 

n (A 1 ,...,A n )= 0f(Ai,..-,A n ), (193) 

e=l 

where 

^ e) (Ai,...,A n ) = T ljl+e -(Ai) 0n-e(A ie - +1 ,...,A in ) 

2<j 2 <-<3e<n 
2<je + l<---<jn<™ 

e 

x ]jT hl (X Jki )g^\X 1} ...,X n ) } (194) 

fci=2 

where gg" (\i, ■ ■ ■ , A n ) represent the coefficients of an arbitrary linear combination. 

(n-l)\ 

At the present stage the number of unknown coefficients for given n and e is r—, rr- 

P & & (n-e)!(e-l)! 

Here we will choose to normalize the n-particle state by the coefficient proportional to the first 
term 7^ )2 (Ai)0 n _i(A 2 , . . . , A„) of the linear combination (I193p . Therefore, the total number of 
coefficients we have to determine is 2 n ~ — 1 which at first sight appears as an impracticable task. 
This number, however, can be truly reduced by invoking the exchange property of the n-particle 
vector under the permutation of the rapidities Ai, • - • , X n . As emphasized in the previous sections 
this is a rather important symmetry present in the explicit construction of the two-particle and 
three-particle states. Inspired by this result we demand that the n-particle ansatz (I193|ll94p should 
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satisfy similar exchange property. For general n the permutation Xj <-> Xj+i for j = 1, . . . , n — 1 
reads, 

^n(Ai, . . . , Aj, Aj + i, . . . , A n ) = 0(Xj, A J+ i)0 n (Ai, . . . , Aj_i, A J+ i, Aj, A J+ 2, • • • , A n ). (195) 

In our analysis of the three-particle vector </>3(Ai, A2, A3) we noticed that the corresponding 
permutation symmetry A2 <-> A3 depends only on the exchange property of the two-particle vector. 
It is expected that this feature extends to general n-particle states providing us a recursive manner 
to solve the problem. From now on we assume that the permutation property (11951) of the n-particle 
vector for j > 2 follows as a consequence of the exchange symmetry of the previous constructed 
(n — l)-particle vector. We stress that this assertion has been explicitly verified up to the four- 
particle state. This implies that relation (11951) should then be valid for each term of the linear 
combination, namely 

0nH^l> • • • ' = ^C\?' • • • 5 Aj_i, Xj+i, Xj, Xj +2 , • • • , A n ), 

for j > 2 and 1 < e < n. (196) 

We now start to explore the consequences of the building blocks permutation symmetry (1 1 96H . 
The first non-trivial case is e = 2 and by substituting <jr n (Ai, . . . , A„) in Eq. (11961) we find, 

g2 +1 \Xi, . . . , X n ) = 9(Xj, Xj + i)g!f\Xi, . . . , Xj-i, X j+ i, Xj, X j+2 , . . . , A n ). (197) 

which can recursively be solved and as result we obtain, 

n 

& y+1) (Ai, . • • , A n ) = ^(A,, A, +1 )^ 2) (A 1; A i+ i, {Xi^j+xl- (198) 

i=2 

In order to make further progress one needs to find the expression for the amplitude 

(2) 

92 (Ai, • • • , A n ). At this point the expressions for the two-particle and three-particle states ob- 
tained in sections 15.21 and 15.31 are of great utility. Direct comparison between Eqs. (I193II194P and 
Eqs.( ll36lll70p for n = 2, 3 reveals us that, 

gf\ Xl , A 2 ) = 1 ^ 2 \X 1 , A 2 ), gi 2) (X 1 , A 2 , A 3 ) = ^^ iJf\Xi, A 2 ). (199) 
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We next note that the i?-matrices prefactor in Eq. (ll99p is directly related to the wanted 
term originated from the commutation rule among the field T^i^) and the (n — 2)-particle 
vector 0„_ 2 (A3, . . . , A n ). Considering this information we are able to infer that the expression for 
di • • • 5 X n ) is, 

gf\xu ■ ■ ■ , x n ) = f[ ^rrw^^ ^ ^ 

and by substituting Eq. (12001) in Eq. (|l98p we finally find, 

#(a 1; ...An) = fl §r^M e <^ A ^)- ( 201 ) 

To emphasize the main points of our procedure we shall also consider explicitly the next case 
e = 3. By substituting 0« (Ai, . . . , A n ) in Eq. (11961) we found that functions g^ 2 (Xi, . . . , A n ) 
satisfy the following constraints, 

93 +1 '^ 3 (Xi, . . . , A n ) = 0(\j, Xj+i)gJ^ 3 \Xi, . . . , Xj-i, Xj+x, Xj, Xj +2 , . . . , A n ), 

for j 2 + 1< j 3 , (202) 

^0'2,j+ ^ _ Xj+ijg^ 2 ' 3 (Xi, . . . , Aj_i, A J+ i, Aj, A J+2 , . . . , A n ), 

for j 2 + 1< j3, (203) 
9 , 3 J '' ?+ ( A i) • • • ) = 0( A i> Xj+xjgJ^ \Xi, . . . , Aj_i, Aj+i, Aj, Aj +2 , . . . , A n ), 

for ./2 • 1 ./:;• (204) 

Once again it is possible to solve Eqs. (l202H204p in a recursive manner and the result is, 

n 

g { t jz) {Xu ...,Xn)= fi ^<( A - X ^)d<(Xi, X j3 )gi 2 ' 3) (X h X j2 ,X j3 , {A, }, , ,). (205) 

1=2 

provided that function #3 3 ^(Ai, . . . , A n ) satisfies the following exchange property, 

# 3) (Ax, . . . , X n ) = 6{X 2 , A 3 )# 3) (Ai, A 3 , A 2 , A 4 , . . . , A n ). (206) 

The remaining task consists in the determination of the fundamental amplitude 
g\ ' (Ai, . . . , A n ) which is done by considering the same kind of arguments employed in the case 
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______ f2 3) 

e = 2. By comparing Eqs. (1193I194|) forn = 3 with Eq. (11701) tell us that function g 3 ' '(Ai, A 2 , A3) = 

(2) 

2T2 (Ai, A2, A3). The prefactors for general values of n are now governed by the first term of the 
commutation rule between the product Tj^i ^2)71,1 (A3) and the creation fields 7l^{Xi) present in 
the leading part of vector n _3(A 4 , . . . , A n ). Collecting all these informations together we find, 

J2,3)/ A \ \ _ TT R(\A2)n R(Xi, A 3 )n \ \\ 

^2)21 ^3)21 



4=4 



We clearly see that Eq. (12071) satisfies the permutation relation (12061) thanks to the exchange 
symmetry of the two-particle function 2^2 A2, A3). The coefficients g$ 2 ■ ■ ■ , A n ) are 

obtained by substituting Eq. (l207p in Eq. (12051) which gives us, 

^'^(Ai, . . . , A n ) = IT , " J 2 sl{ x 3 ' 3 N2^<(^i) Ay 2 )^<(Ai, Aj 3 ) 2 ^ r 2 (2) (^l> ^2) ^ja)- ( 208 ) 

7=2 K V A i> A j2>21 A h)21 

Let us now generalize the above results for arbitrary index e. The analysis of the consequences 
of the permutation symmetry (I196p is made by using mathematical induction, leading us to the 
following general constraint, 



# ,/: (A; a„) = n ]jo ,: ( A '-./2 ./.•{A;}:,t:':.,,)- ( 209 ) 

i=2 1=2 



We are again left to determine the expression of the basic function gf' '"' e ^(Ai, . . . , A n ). To 
make progress in this direction we first pause to discuss the results we have so far obtained. From 
Eqs. (l200ll207p we note that functions g^ (Ai, . . . , A n ) and ^'^(Ai, . . . , A n ) are given in terms of 
product among universal prefactors with the off-shell amplitudes \T\ (Ai, A2) and iT\ (Ai, A2, A3). 
The former amplitude is a central element of the eigenvalue problem for the one-particle state 
while the latter play analogous role for the two-particle state. It is natural to assume that the 
structure of functions gf' '"' e ^(Ai, . . . , A n ) for arbitrary e will follow the same pattern noted above. 
In other words, that these elementary functions can be fixed after the solution of the eigen- 
value problem for the previous (e — l)-particle state is completed. In analogy to sections I5.1H5.3I 
such solution will provides the expression for the respective amplitude g-i^i^^Ai, Xj 2 , . . . , A Jg ). 
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The expression for g^'' ■ ■ ■ , A n ) is then obtained by multiplying prefactors involving the ra- 

tios R(\i, Xj)\'\/ R(Xi 7 Xj)l'\ with the amplitude e-i^g-i(Xi, Aj 2 , . . . , A Jg ). As before the prefactors 
are easily accounted through the commutation of the product Ti^X?) • • . 7i 5 i(Ae) with operator 
n _e(A g +i, . . . , A n ). Considering this reasoning we obtain, 

gt-'^Xu , K) = f[ fl ^TW^i^ ■ ■ ■ > ( 21 °) 

i=e+lj=2 V *' J/21 



and by substituting Eq. (12101) in Eq. (j209p we found, 

^'•••' ie) (Al, ...,A n ) = TT Y\ jrpr J k \ 2 l #<(Afc 2 , Aj fcl )e-l^g?i(Al, Aj 2 , . . . , Aj g ). (211) 

fe 1= 2 fc 2 =2 ^^^21 

We have reduced the determination of an arbitrary n-particle vector to the computation of 
only N — 1 coefficients e-i^-i(Xx, A J2 , . . . , Xj s ) for e = 2, . . . , iV. Collecting together Eqs. fll93[ 
11941) and Eq. (1211j) we are able to propose that an educated ansatz for the n-particle vector should 
be, 

m(n,N-l) 

^l,l+e(Al) 2^ n _g(Aj- +1 , . . . , AjJ e-lf =-l(Al, Xj 2 , . . . , AjJ 

g=l 2< J2 <-<ie<n 

2 <ie+l<---<jn<n 

e n _^(_\ . _\ . "j 11 

x n r ^(Ai fel ) n Rix^ x P <{x ^ x ^ ] (2i2) 

fcl=2 k 2 =e+l ^ Jfc 2' 3fci/21 

where the normalization qJq(X) = 1 is assumed. 

Having at hand a proposal for the n-particle vector we can turn our attention to the action 
of the diagonal operators T a)a (X) on the multi-particle state |$ n ) = n (Ai, . . . , A n ) |0). This study 
is fundamental to find the expressions of the off-shell functions e-i^g_i(Ai, Xj 2 , . . . , A Jg ) in terms 
the i?-matrix elements. It should be stressed that the approach put forward here is clearly self- 
consistent. After the solution of the eigenvalue problem for the sector with (n — 1) particles 
we have the basic ingredients to pursue the next sector having one more particle and so forth. 
Considering our previous results for one-particle, two-particle and three-particle together with the 
help of mathematical induction we are able to determine the action of the diagonal fields on |$ n ). 
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The final results are, 



T a>a {\) |$ n ) 



X 



n m{a— l,t] 

.(A) n ft(A, A t ) i*„> - x; e r a- Pl « +t - P (A) 

i=l t=l p=M{0,a+t-N} 

i=l,...,n \ sr-(a-p) 



£ 



i<ji<---<j(i_ p) <™ 

l<3(t-p+i)<-<Jt<n 
t-p 



fe=l 



0n-t({Aj} i ^{' ) '_'_'_"- t ) t-pJF t (A, Aj 15 . . . , Aj t ) 

\ 
/ 



i=1 -R(Ai, Aj fc ) 2 |i 



j=t-p+i 



/ 



n n <( A iP A ^)i°) 
/c=i ;=t-p+i 



(213) 



The structure of the recurrence relations for the off-shell amplitudes C .F 6 ^(A, Ai, . . . , A&) is as 
follows. For c 7^ and c ^ b our findings for the two-particle and three-particle states given in 
Eqs. (ll44)181jll82p are sufficient to guide us to the following general structure, 

b c 



R(\i, i 



c T^ a \\ Ai, . . . , A fe ) - o^ r 6 ( " c (A, A( c+ i), . . . , Xb)cFc a+b Ai, . . . , A c ) ] [ | f ( •> , 

i=c+lj=l K ( A i' A j >2,1 



for 6 = 2,...,iV-l: a = 1, . . . , N - b: c = 1, . . . , b - 1. 



(214) 



The expressions for c = 0, b are, however, more complicated and their multi-particle extension 
required us to perform explicit computations up to the four-particle state. Considering the help 
of mathematical induction as well as the expected properties of these functions under exchange of 
their rapidities, see for instance Eqs. (jl46[ll84p . we are able to infer that, 



o^F b a \x, Ai, . . . , Aft) 



£ 



R(X,X 



lJa+e,l 



1 -^(A, 



E 



O-T^g^A, Aji, . . . , A(ft_e)) 



2<J 1 <-<i( b _g)<!> 



b-e 6-1 



X e-l-? 7 e-l(Al) A 



J(6-e+l) ' " " " ' 



n~i — »- R(^\j^ , A^ ) | | 



i 1= l ; 2 =6-g+i ^( A Ji! ' A j; 2 ^1 
for 6 = l,...,iV-l; a = 1, . . . , N - b 



(215) 
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and 



6 ^ a (A, Ai, 



6-1 



- ~E E f^b a 

f=0 l<h<h<-<l<j,-j)<b 

nTT a t\ \ Rftl.i ^)|i 

^ i\ e<(A -' A '- ) fl(A„A I ,)I;: fl (A I „A,)!;: 

for 6 = 1 JV-1: a = l N — b. 



(b-f) 



(216) 



Here we stress that the recurrence relations (I214II215|216I) have as initial condition the overall 
normalization o-^o (A) = 1 as well as the expressions for the one-particle off-shell amplitudes 
o^ B) (A,^) and !^ o) (A,/i) defined by Eq.JT2SJ. 

The solution of the eigenvalue problem for the multiparticle state is completed by summing 



over Eq. (12131) for a = 1, . . . 

b Tr' (A, Ai, . . . , A 6 ) we find, 



N. By substituting the expression (I216p for the off-shell amplitudes 



N 



n N-t 



t-1 



t(x) i*„> = w ^ II p ^ ^ i*»> -EE T ^( A ) E E 



a=l 



8=1 



t=l a=l 



p=0 i<ji <•■•<«><« 
l<i( P +l) <-<it<n 



(a). 



V 

t p n 

n n n 

r=p+l s=l 



i=1 -R(Ai, A Js ) 2 ' 1 



J 



i=l 
i&U—jt 



/ L 



r=p+l 



<=1 ^(Aj r ) Ai) 2 'i 



X 



0Wr,*3.)- 11 ^i( A >) 11 pn : x2i 11 



s=l 



r=p+l 



i=1 Aj r ) 2 |i s= i R(^j s , -V^i ^(A>, AjJi^i 

i&u—jt 



• (217) 



The unwanted terms are canceled out by requiring that functions inside the bracket of Eq. (12171) 
are null for t — 1, . . . , m(n, N — 1). The situation here is exactly the same found already for the 
two-particle and three-particle state. It is sufficient to consider the constraint coming from the 
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case t = 1 which leads us to the following Bethe ansatz equations 

il.l r>f\ \ \2,1 



-^#4 - f[ d ( X h> X i) RiKAi) ii ^ Xh) i\ for Jl = l ; ...,n. (218) 



i=l 



To verify the cancellations for the remaining values of index t one just has to substitute Eq. (12181) 
on the bracket of Eq. (l217p and to consider the exchange property (I140p of functions 8(Xi, Xj). The 
elimination of all unwanted terms by Eq. (12181) implies that the n-particle eigenstate is, 

N n 

A n (A) = MX) J] p -( x > x i)- ( 219 ) 

a=l i=l 

We close by mentioning that all formulae of this section pass by the following tests. For 
iV = 2 we recover the well known algebraic Bethe ansatz solution of the six- vertex model |10j . 
In this special case one has to note that both numerator and denominator of Eq. (11381) vanish 
and the respective limit gives us #(A,/i) = R(X, /x)^/ R(X, /-Oi'i- For iV = 3 we reproduce the 
algebraic Bethe ansatz construction proposed by Tarasov [12] for the nineteen-vertex models. 
Strictly speaking, our results should be seen as extensions of the above mentioned works since 
no assumption of spectral parameter dependence for the i?-matrix has been made. We also 
verified that the expressions (I218|219 l ) reproduce the Bethe ansatz equations and eigenvalues of 
the higher spin generalization of the six- vertex model [IT] . We have also computed explicitly 
the most complicated off-shell amplitudes of the four-particle state such as o^i (A, Ai, . . . , A4) 
and 4^4°^ (A, Ax, . . . , A4). These computations have been checked to be in accordance with the 
recurrence relations (I215ll216p . 



6 Conclusion 

In this paper we have shown how the algebraic Bethe ansatz method works for arbitrary vertex 
models whose i?-matrix commutes with one U(l) symmetry. This invariance guarantees the 
existence of a reference state in which the respective monodromy matrix acts triangularly. We 
recall that this is necessary requirement to start performing an algebraic Bethe ansatz analysis. 
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The transfer matrix eigenvalue problem is viewed under the general perspective that its solution 
should not depend on a specific functional form of the i?-matrix. We argued that the algebraic 
formulation of the Bethe states of the transfer matrix can be done solely on basis of the Yang- 
Baxter algebra, the Yang-Baxter equation and the unitarity property satisfied by the i?-matrix. 
In fact, we have presented a method to obtain the appropriate commutation rules between the 
monodromy matrix elements for arbitrary number N of edge states. Moreover, the necessary 
identities among the i?-matrix elements to solve the eigenvalue problem are derived from the 
respective Yang-Baxter and unitary relations. 

This approach provided us the expressions for the on-shell properties such as the transfer matrix 
eigenvalues and the Bethe ansatz equations as well as the structure of the off-shell Bethe vectors 
in terms of the arbitrary i?-matrix weights. The respective off-shell amplitudes are determined 
by means of a recurrence formula whose input are special entries of the i?-matrix. Note that the 
previous understanding of off-shell structure of U(l) integrable models was restricted to the six- 
vertex pU] and nineteen- vertex models [12J. By contrast, this paper offers us the basic ingredients 
to compute such properties for any U(l) vertex model. In a forthcoming paper [20] we shall indeed 
exhibit the explicit expressions of the off-shell amplitudes of some U(l) invariant vertex models 
such as the classic higher spin XX Z — S Heisenberg system and the non-compact vertex model 
based on the SL(2, R) algebra. Interesting enough, we found that all the off-shell amplitudes 
C .F 6 (A, Ai, . . . , Aft) factorize in terms of product of elementary functions. 

We hope that the framework developed in this paper will also be relevant to solve other 
families of integrable models by the algebraic Bethe ansatz. In particular, that the universal 
structure of our expressions for the eigenvectors and the off-shell amplitudes should play the role 
of cornerstones in nested Bethe ansatz solutions of integrable models with i?-matrices commuting 
with more than one U(l) symmetry. This idea has indeed worked for the Hubbard model whose 
algebraic Bethe ansatz solution [19] encoded the main features found in the algebraic Bethe ansatz 
formulation of the N — 3 U(l) vertex model. In general, we expect that the algebraic solution of 
a iV-state integrable model with m distinct U(l) conserved charges will be described in terms of 
a (N — m+ 1) x (N — m + 1) monodromy matrix. The partition of the monodromy matrix from 
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N x N to (N — m + 1) x (N — m+1) are going to be dictated by the form of the U(l) operators that 
commute with the i?-matrix. This means that the respective length of the recurrence relations for 
the eigenvectors and off-shell amplitudes will be governed by the effective number N — N — m + 1. 
Now the off-shell amplitudes C T\ , (A, Ai,...,A&) will be seen as vectors with m 2 components 
while 8(X,fi) will behave as a m x m auxiliary factorized R(X, fi) .R-matrix satisfying the unitary 
property. We recall that this scenario is known to work for specific vertex models with N — 2 
U(l) symmetries such as those based on the non-exceptional Lie algebras [13J. This has been also 
verified for the vertex model related to the fundamental representation of the U q [G2\ algebra |21j . 
This system is rather involved since the size of its R-matrix is 7 x 7 while its algebraic Bethe ansatz 
solution is formulated by a 5 x 5 monodromy matrix. In all these cases the corresponding auxiliary 
^-matrices -R(A, /j,) is derived from the analysis of the commutation relations following the general 
method explained in section [31 Their factorization and that they fulfill the unitarity relation 
were verified using the respective specific weights. The results of this paper suggest that such 
properties can be derived as a consequence of the Yang-Baxter equation and the unitary relation 
of the original i?-matrix we have started with. Hopefully, this could be shown by generalizing the 
method discussed in section HI 

Finally, we think that this paper highlights the expectation that the Bethe ansatz properties 
of integrable models should be exhibited in terms of universal formulae depending only on the 
i?-matrix amplitudes. We hope that our results will inspire the search of a general recipe to 
solve integrable systems invariant by many U(l) symmetries through the algebraic Bethe ansatz 
method. 

Appendix A: The two-particle state 

In this appendix we provide the technical details entering the solution of the two-particle eigenvalue 
problem. Here we shall assume that the easy unwanted terms mentioned in section \572\ have been 
already canceled out. Therefore, the results to be described in what follows concern with the 
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symmetric two-particle ansatz, 



1*2) = Tj. 2 (A, )7]. 2 (Ao) |0) - ^^g«;i(A 2 )T 1 ,3(A 1 ) |0) . 



(A.l) 



#(Al, ^2)3^1 

The elimination of the easy unwanted terms brings considerable simplifications on the action 
of the diagonal operators T a;a (X) on |$ 2 ). Considering all the steps mentioned in section ISTiZl we 
find that T a>a (X) |$ 2 ) has to be divided in five different parts. We find that they are given by, 



T a , a (A)|$ 2 ) = w a (X)P a (X,X 1 ) P a (A,A 2 ) ^ 1) ) + P a (A,A 1 ,A 2 ) 



(2) 



8a E ^i(^)i^ a) (A, Ai, A 2 |i)T a)a+1 (A)T 1)2 (A i ) |0) 



i,j=l 



51 w 2 (X l ) H[ a - 1) (X, A 1; A 2 |i)T a _ lia (A)T li2 (A J ) |0) 



i.3=l 



^- 1 ' N w 1 (X 1 )w 1 (X 2 ) 2 ^ a) (X,Xi,X2)T a , a+2 (X) |0) 
2 

Sl' N E ^i(Ai)^ 2 (A i )i^°~ 1) (A, A x , A 2 |z)T a _ 1>a+1 (A) |0) 



^ 2 (Ai)«; 2 (A 2 )o^ 2) (A, A 1; A 2 )T Q _ 2 , a (A) |0) 

for 1 < a < N, 



(A.2) 



Here we recall that functions P a (A,/i) have been defined in the one-particle problem, see 
Eq. fll27p . The dependence of the extra functions P a (A, Ai, A 2 ) on the P-matrix elements are given 
by, 



Pi(A,A l5 A 2 ) 



R(Xi, A 2 )|} 



R(Xi, X 2 )lf 



R(^2,X)l' 2 R(Xx,X)ti R(Xx, X2)s'lR(Xi, X) 2 \ 



\2,2 



2.2 



\2,1 



-R(A 2 , A) 2 'i R(Xi, X)l'[ R(Xi, X 2 )'l'[ R(Xi, X)l'[ 



3,1 



\3,1 



(A.3) 



Pi(A, Ai, A 2 ) 



R(Xl, A 2 )g^ 

-R(Ai, X 2 )l'l 



r(\ ^y a l hl Dt 1 ^ 



(A,A 1 ) + P(A 1 ,A 2 )^^(A,A 1 ) 



_ R(K A 2 )SS /Jf +1 ' 0) (A, Ai) R(X 1 , X 2 f 3 j D^°\x, Ai 

^(A.AQ^A,)- 1 ' ^ 

— r-?v^ TTi - i° r 2 < a < N — 2 

^(A,Ai) i2(A,A a );} 



(A.4) 



75 



Pn-i(X, Ai, A 2 ) 



-R(Ai, X 2 )^) 



2.2 



-R(Ai, A 2 , 3 :i 
i 

i?(A,A 



R(\ A 2 )^ 1 1 ' 2 R(X, Ai)^ 13 i?(Ai, A 2 )^ -R(A, Ai; A | 



vAT,2 



AT,1 



-R(A, A 2 )jv'i 



-R(A, Ai)jv2 



-R(Ai, A 2 ) 3 \ i?(A, Ai 



N,2 



X 



-R(A, Ai)^, 1 ' 3 



^(A, ^l)jv!2 



UjV-1,3 it l A ' A lJiV-l,3 



-R(A, Ai)^! 1 ' 2 



-R(A, Ai)^'J 



-R(A, Ai)^_}'2 -R(A, Ai)^ 12 



.D^" 1 ' 1 ^, Ai) R(X, X 2 )nZ 2 1 ' 2 1 
Di N - 1 >°\x,X 1 )R(X,X 2 f N Z{j 



(A.5) 



Pn(X, Ai, A 2 , 



-R(Ai, A 2 ) 3 | 1 

-R(Al, A 2 )g'l 



A 2 ) 3 2 -R(A, Ai 



N,3 



R(X, Ai)^' 2 ! 3 R(X, X 



2)N,1 



R(Xi, A 2 ) 3 ^ i?(A, Ai)^' 2 



-^(A, ^i)tv!2 



-R(A, ^2)jv 1 



(A.6) 



For sake of clarity we have divided the functions proportional to the unwanted terms in two 
distinct categories. The first class consists of those that carry an internal discrete index dependence 



cn { b a \X, Ai, A 2 |l) or c H ( b a \X, Ai, A 2 |2) where 1 < b < 2, 1 < a < N - b, b - 1 < c < 1. The 
expressions for functions c Ti b (X, Ai, A 2 |l) are 

1 wf ) (A,Ai,A 2 |l) 



/(a) 



i?(A 2 ,Ai)i; i --{a),, \ \ r , 

2\\-Fi (A, AiJ for a = l, 

-R(A 2 , Ai) 2 ' 1 



AT 



wS a) (A,Ai,A 2 |l) = P 2 (Ai,A 2 ) ^r(A,A 1 ) for a = 1, . . . , N - 1 



(«), 



for a = 2, 



AT 



(A.7) 

(A.8) 
(A.9) 



iW^&Ai.Aall) = i^CA.AOoM^CA.Ai.Aa^) 

The explicit expressions for functions c TC b (X, Ai, A 2 |2) are in general very cumbersome. For- 
tunately, thanks to the exchange symmetry of the two-particle vector (11391) there exists a direct 
relationship between functions c H b a \X, Ai, A 2 |2) and c w[' 1 (A, Ai, A 2 |l). In fact, by considering the 
permutation Ai <-> A 2 on the two-particle state result (1A.2|) we are able to derive the following 
consistency relation, 

c H { b a) (X, A 1; A 2 |2) = 0(A 1; A 2 )X a) (A, A 2 , A^l) (A.10) 
However, in the course of our analysis we find that explicit expressions for the following special 



families iH[ a) (A, Ai, A 2 |2) and Hf^ 1 (X,X 1 ,X 2 \2) for a = 1, • • • , N - 1 are indeed very relevant. 



/(<*+!) 



For this reason it is necessary to quote their expressions, namely 



n?(x,x u x 2 \2) 



R{X 2 , A) 2 ' 2 



-R(Ai, A) 2 | 2 i?(A 2 , Xi)f 2 R(Xi, X 2 )i]i R(Xi, X) 2 



v2,l 



2.2 



-R(A2, A) 2 j 



^ 2 (Ai,A) 



\3,1 



-R(Ai, A) a|i -R(A 2 , Ai)|} -R(Ai, A 2 ) 3 ^ i?(Ai, A) 3 ;} 

(A.ll) 
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1 H! a) (A,A 1 ,A 2 |2) 



Dj afi \\,\i) fl(A,A 2 )^ M R(\,\i)£ _ fl(A,A 1 )&,ifl(A2,A 1 )S 
D^ +m (A, Ai) i2(A, A 2 )S; : } i?(A, A^i i?(A, AOS};} i?(A 2 , AO?;} 



^(Ai,A 2 )^ 



-R(A, Ai)°+ 2j i 
-R(A, Ai)^2a 



-^(A, Ai)„+ lj2 
-R(A, Ai)^2 



-R(Ai, A 2 ) 3 ;i 



^(A,Ai)°+2^ ^(A,Ai)°+^2 
R(x,Xi)lXti R(K Ai)"+^2 



for a = 2,..., N- 2 (A.12) 



oM a) (A,Ai,A 2 |2) = 



R(X, A 2 )o+i j i 



-R(A, Ai)"'i 



-R(A, \i)l+n R(Xi, A 2 )oi 



-R(A, A 2 )°+l^ -R(A, Ai)"+J|i i?(A, Ai)^^| -R(Ai, A 2 )|} 



(A.13) 



The remaining functions to be defined are o^^^A, Ai, A 2 ) and 2 jF 2 a - ) (A, Ai, A 2 ). Their explicit 
expression in terms of the weights are, 

-R(A, Ai)"^ 1;1 R(X, A 2 )^i -R(A, Ai)„+ 2 ,i -R(Ai, A 2 )^ 



•7~2 ^(A, Ai, A 2 ) — 



R(X, Xi)lXti ^(A> A2)a+^i -R(A, Ai)^+2^ -R(Ai, A 2 )l'i 

for 1< a < N - 2 



(A.14) 



2^ 7 2 1 ' ) (A, Ai, A 2 ) 



(A, Ai, A 2 ) 



fi(A 2 , A) i 
i?(A 2 ,A) 



2,1 



£>2 2,1) (Ai,A) 



^(Ai, A)^3 i?(Ai, A 2 )g 



2,2 



R(Xi, X)l'\ R(Xi, A 2 ) 



2,1 



3,1 
3,1 



-R(Ai, A 2 )ai 



-R(A, Ai)^ 2 ,i 
-R(A, Ai)^+2'i 



-R(A, Ai)a+ 1;2 
R(X, Xi)lXi% 



R(Xl, A 2 )g'i 



-R(A, Ai)^i ^(A, Ai)"^ 



-^(A, Ai)"+ 2 'i 



-R(A, Ai)£li'2 



-R(A, A 2 )o+i i 



-^(A, Ai)^'i 



-R(A, Ai)°+ 2 ^ 



^ +1 ' 0) (A, Ai) R(X, A 2 )SJ;J i?(A, A0SJ;} i?(A, Ax)^ 

for 2 < a < N - 2 



(A.15) 



(A.16) 



To make progress on the structure of the two-particle problem we first investigate the wanted 
terms. It is essential to assure that the functional form of the products proportional to either 



(i) 



or 



$ 2 2 ^ ) are exactly the same. This means that we have to demonstrate the following 
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relation, 

P a (A,A l! A 2 ) = P a (A,A 2 ) for a = l,...,N. (A.17) 

The factorization property ( 1A.17I) follows from a series of identities we have derived in section 
(14. 2p . The particular cases a = 1 and a = N are an immediate consequence of the Yang-Baxter 
equation. This is easily seen by comparing the definitions of the corresponding P a (X, Ai, A2) given 
in Eqs.( lA.3|A.6[) with the identities ( 19611971) . respectively. In order to show property ( 1A.17I) for the 
remaining cases 2 < a < N — lwe have to combine together the identities derived from both the 
unitarity relation and the Yang-Baxter equation. In order to show Eq. (lA.17l) for 2 < a < N — 2 
we just have to replace the determinants D±\\, Ai) and D§' 2 \\, Ai) in the Yang-Baxter identity 
(11141) by using the relations ( 1931I9T1) derived from the unitarity property. The proof for the case 
a = N — 1 requires to consider the extension of Eqs. (l93f9ip when the index i goes to the value 
N — 1. As a result of such analytical continuation we found, 



lim 



D 



(A, Ai 



TV-1 £)0+1.0) 



(A, Ai 



(A, Ai) 



and 



jff 0) (A,AQ 
If - l D^ ai (X,X 1 ) 



lim 



-R(A, Ai)^^ 1 ' 3 



^(A, Ai)^'2 



TV,2 
UjV-1,3 



R{\,\ 



(A, Ai 



-R(A, ^i)na 



(A.18) 



R(X, Ai)^ 1 ' 2 



R(\ Ai)^'i 
R{\ Ai)^ 1>2 



R(\Xi) 



TV, 2 
TV, 2 



(TV+1,4) 



(A, Ax) 



Dl N > 3) (\,\ 1 )D. 



(JV+1,3) 



(A, Ai) 



(A.19) 



To complete the proof we have to consider the analytical continuation of Eq. ( fTT4l for a = N — 1 
by means of relations (I115|lll6p . The replacement of the determinants D\ (A, Ai) and D^' 2 \\, Ai) 
in this analytical extension of Eq. (11141) with the help of Eqs. (I93IIA. 18|IA.19I) leads us to relation 
( !A~T7l) for a = N - 1. 

We now turn our attention to the structure of the unwanted terms. In order to obtain the result 
(11411) presented in section l5\2l we first identify functions iH^\\, Ai, A2 1 1) and (A, Ai, A 2 ), 



i7i 2 (A, Ai, A 2 |l) — \T\ ''(A, Ai, A 2 



(A.20) 
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The next step is to carry out some simplifications on functions C T 2 (X, Ai, A 2 ) for c = 0, 1, 2 
proportional to undesirable terms with spin s — 2. It is possible to express these functions in 
terms of recurrence relations involving the one-particle weights o^i^A, Ai, A2), i-F^^A, Ai, A 



2 



and certain i?-matrix amplitudes. For function 0^2 (X, Ai, A2) this follows directly by using the 
definition (11281) in Eq. (1A.14|) . namely 



R(X, ^l)o+l,l T (a+1) M x \ -R(^) ^l)a+2,l 

— — — ^+2jo>i (a, a 2 j + 

^(A, AlJa+2,1 K { A , A l)a+2,l 



J- 2 [A, Al, A 2 J — X a+2,1°"M (A, A 2 J + xa+2,ll^ 1 l A l' A 2j- l A .Zl y 



In the case of 1J2 (A, Ai, A2) we have to use the identification (IA.20I) together with the defi- 
nitions (1A.8|A.9j) as well as the property (lA.lOj) . After few manipulations we find 



1 ^ ) (A,A 1 ,A 2 ) =o^ ) (A,A 2 ) 1 ^r i) (A,A0^4 1 ^- (A.22) 

where we have used the relation -R(A 2 , Ai)2^P 2 (A2, Ai)0(Ai, A 2 ) = -R(A 2 , Ai)^} which can be derived 
considering Eqs. fll27lll38lfT40]) . 

The last simplification concerns to show that the definitions flA.151 IA.16j) for functions 
2 J 7 2 a ' ) (A, Ai, A 2 ) are equivalent to the following expression, 

2 R(X ■ A-)^'} R(X A ) 1 '^ -1 

2F 2 (X, Ai, A 2 ) = -qF 2 (A, Ai, A 2 ) - ^2 x ^2 (a, Ai, Xj) 1 ' — — — 9(X k , Xj). 

i,j=i K \ A j' AiJi t i K[Ai, Aj) 2jl k=1 

(A.23) 

In order to show the equivalence between Eqs. (lA.15|A.16p and 0A.230 it is necessary to use 
particular relations coming from the symmetry relation Eq.f lA.10j) . namely 

oM a+1) (A,Ai,A 2 |2) =e(Ai,A 2 )oM a+1) (A,A 2 ,Ai|l) for 1 < a < N - 2 (A.24) 



iHf j (A,Ai,A 2 |2) = 0(Ai,A 2 )i^r(A,A 2 ,Ai|l) for 1 < a < N - 2. (A.25) 

We first substitute the expressions for these functions given by Eqs. (1A.7IIA.1 If A. 12[) in the 
relation (IA.250 . We then reorder the corresponding amplitude ratios of type pl^'^i for (A, jj) = 
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(A2, Ai) and (A,//) = (Ai, A) with the help of identity (179|) . As a result we obtain, 

R{X2, X)i2 p ^ _ -^O^i R(Xi, AQ^i _ #(Ai> Ap^i -R(Ai, A)^ 
i?(A 2 ,A)2;} ii!(A, Ai)a|i R(Xi, A 2 )|i i?(Ai, A 2 )|} -R(Ai, A)|} 

= fl(A a , X 2 ) ^ U ^\ 1 Ti 1 \X, A 2 ) for a = 1 
-cl(,Ai, A 2 J 2 'i 

(A.26) 



(0,0) 



(a+1,0) 



(A, Ai) _(a)/, , x -R(A, A 

l)a,l -r( a ) / \ \ 

R{Xi, A 2 ) 2 'i 

(v 1 (A A 2j — , , a+1 x - o-/~ 1 (A, AlJ 



(A, Ai 



i?(A,A 



-R(Ai, A 2 ) 3 ^ 



-R(A, Ai)"+ 21 



^l)a+l,2 

R(X, Ai)^i' 2 



-R(Ai, A 2 )l'i 



R(X, Ai)^ 2 ^ 

r(x, Xi)IXI'a 



R(X, Ai)^2 
R(X, Ai)^ 2 



9(X 1 ,X, 



R(X±, X 2 )l , \ 



R( X ^ X 2)lj (a) 

•".(Al, A 2 J 2 ,i 



for 2 < a < A - 2. 



(A.27) 



We next multiply Eq. flA~27l) by function ^i (a+1) (A, Xi) and Eq. (lA~24l) by function o^i (a) (A, Ai). 
By subtracting the former equation from the latter equation and by considering the explicit 
expressions for H ( f' +1 \X, Ai, A 2 |2) given in Eq. (lA.13p we found, 



yr\ (A, Aijo^i (A A 2 J +2 ,i 



0^1^ (A, A2)o-? r i a+i ' , (A, Ai) 



(a+1,0) 



-R(A, Ai)"+}^ 



x ^ (a+1) (A,A 



-R(Ai, A 2 )||i 
i2(Ai, A 



34 (H 1 {A, Ai J 



(A, A x ) 

R(X, Ai)^ 2 ,i ^(A,Ai)*^ 1(2 
-R(A, Ai)^+oi ^(A,Ai)„ti' 2 



2/3,1 



i?(A,A0: 



a+1,2 0/ \ \ \a+l,2 
-H-IA AiJ a+1 2 



+2,1 

+2,1 
+2,1 



-R(A, Ai)*+i' 2 



R(X2, Ai)^ a 
-R(A2, Ai)^} 



qT[ (A, Ai) 



-n-(Ax, A 2 j2,i 



(A.28) 
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After few manipulations in Eg. (1A.28I) we are able to write it as, 



(A, 



D 



(a+1,0) 



-o^i (A, A 2 )o^ r i + (A, Ai' 



R(\ Ai)°^ 
-R(A, Ai)o +1 'i 



-^(A, Ai)"+ 2 ,i ^(A, Ai)„+ 1)2 
i2(A,A 



\a+l,2 td(\ \ \a+l,2 
lJo+2,1 K V A > A lJa+l,2 



•R(Ai,A 2 ) 3 |? #(Ai, A 2 )|?i2(A, Ai)a 



+2,1 



-R(Ai, A 2 )|} -R(Ai, A 2 )g'i -R(A, Ai)^ 2 ^ 



3,1 



\a+2,l 



(a+l) 



-R(A, Ai)^i ^(A, Ai)a+{| 2 

-^(A, Al)„+ 2 'l R(\ Al)a+i^ 2 
-.(a). 



l A ' A 2J , Na +2,1 ~~ 0^ 1 (A A ly 



-R(A 2 , Ai)^} 

-R(A 2 , Ai) 2 'l 



X!^ a+1) (A,A 2 ) 



fi ( Al ' A2 )M 77(a) x , -r 
— — 2j ^ 1 l A , A 2jl^" 1 

-fH A l, A 2j 2 1 



0(\i,\ 2 ) 



R(X, Ai; a+21 

{a+1) (X,Xi). 



(A.29) 



which in fact shows that Eqs. (lA.16IIA.23p are equivalent once we consider Eq. (lA.22p . 

Finally, it remains to show that Eg. (IA.23I) for a = 1 is equivalent to Eg. (lA.15j) . In this case, we 
multiply Eq. flA~24j) for a=l by o-^f^A, Ai) and Eq. flA~26|) by o^f 5 (A, Ai). By subtracting these 
equations we found, 

-R(A, Ai) 2 'J R(X 2 ,X)li 



o«^i ''(A, A 2 )o^i ''(A, Ai) 
+ o-^j ^(Ai, A 2 )o^ r i ^(A, A 



jP 2 (Ai, A)o^ 2 ' ) (A, Ai 
i 



-^(A, Ai)|} i?(A 2 ,A) 2 j 



2 o i2(A 2 ,A 1 ) 1 ; 1 (2) (i) N 

jo-' i ( A ; M)o-^i ( A , A i, 



^(AijA),! i?(A 2 ,Ai) 2 '| 



#(Ai, A 2 )^t; — - ryi i-^j^CA, A 2 )o^ 2 ^(A, Ai). 



i2(Ai, A 



(A.30) 



2^2,1 



We have now to eliminate the terms P 2 (X U A) and ^i (2) (A, Ai)§^|t from Eq. (fA~30l) . This 
step is implemented with the help of two identities derived from the unitarity relation. More 
specifically, these are the identity (185]) and that coming from Z/[l, 3]f of Eq. (!78|) . namely 



R(\, fi^Rifi, A) ji + R(X, i*)%R(ii, A)$ + R(X, riftRfa A) ft = 0. 



\3,1 



\2,2 



3,1 



\3,1 



\3,1 



(A.31) 
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and as result we obtain, 



#(A 2 ,A)?; 2 (2,i), x iU x , R ( A ^ A )lj _ fl(A 2 , AQ^ (2) (1) ^ 



-R(A 2 ,A) 2 \ i?(Ai,A) 3 ' 1 i?(A 2 ,Ai) 2 ' 1 



#(Ai,A 2 )—— — .2\iF 1 (A, A 2 ) ^ r 1 (A, Ai)— o^ 7 ! (A, A 2 ) -3^ 
-KlAi, a 2 j 2 ; 1 K(A, AiJ 3)1 



^(A, A x ) 3 ;i 



+ o^i (2) (Ai,A 2 )^^g. (A.32) 



which shows the equivalence of Eqs. (lA.15IIA.23j) for a = 1. 

Putting all these results together we will find the final results (I14ip presented in section (15. 2p . 



Appendix B: The three-particle state 

In this appendix we describe some of the technical details used to solve the three-particle problem. 
The steps necessary for the last two terms of the three-particle vector (11701) have already been 
discussed in the text of section 15.31 Therefore we shall concentrate our attention on the most 
intricate part of the calculations which is associated to the first term / 7i ;2 (Ai)0 2 (A 2 , A3). In this 



case we have to carry the operators ^^(A) through the state 



(i) 



T lj2 (Ai)0 2 (A 2 ,A 3 )iO>. The 



first step is done with the help of Eqs. ( 12711291) and as result we obtain, 



21.1(A) 



(i) 



R(X 1 , A)^ 



-R(Ai, a)1'1 



-R(Ai, A) 2 ^ 



^MAOTx^A) 



R(\i, A) 2 \ 



tTi, 2 (A)T 1)1 (A 1 



2 (A 2 ,A 3 )|O>, 



(B.l) 



(1) 



a+1 



^ 0) (A,A 1 )T li2 (A 1 )T a , a (A) + ^^ e "- 2) (A,A 1 )T 1 , e -(A 1 )r a , a+2 _ e -(A) 



a Pf\ \ \ a ' 2 PA \ \e,a-e+l 
K ^ A ' A l )a+l,l K \ A > A l)a,l 



i)o+i!i ^(A> Ai)a'i 



o,l 2e, a +l 



o-l 



- £ 



e=l 



-^(A, Ai) a \ 



-^l ^g,a(A)7^_g +12 (Ai y 



e=3 

(A)T a _ g+1 , 1 (A 1 ) 
2 (A 2 ,A 3 ) |0), for 2<a<iV-l, (B.2) 
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T N , N (X) <5> ( 3 



) 




-R(A, Ai) JV ' +2 _g i g 
-R(A, Ai)^'} 



^1 ,g ( A i ) 7/v,7V+2-e ( A) 




—I R(X, Al)^ 

^ i2(A,A!) 




7;- jV (A)T Ar _ e - +li2 (A 1 ) 2 (A 2 , A 3 ) |0) . 



(B.3) 



By examining Eqs. (1B.lHB.3l) we conclude that our next problem is to compute the action of 
operators T d+a a (X) with d = 0, . . . , N — a on the two-particle state </> 2 (A 2 , A 3 ) |0). The main stages 
to disentangle this problem is as follows. For the diagonal fields d = this task is implemented 
by using the results ( 1141ft for the two-particle state. This operation, however, is able to pro- 
duce products of creation operators of the form 7i ) 2(Ai)7^ 1)0l+ i(A) / 7i )2 (A J ) |0) with j = 2, 3 and 
^i,2(Ai)7^ 2ia2+2 (A) |0) for cii = a — i, . . . , a that need to be reordered as far as the rapidity A is 
concerned. The product of the creation fields is sorted out by commuting the operators 7i i2 (Ai) 
and T ai!ai+ i(\) by using Eqs.(j46|47j). After that we carry on the diagonal and annihilation opera- 
tors through the field Ti 2 {\j) with the help of Eq. fj52ll26j) . These two steps together are able to 
generate the following product of creation fields 7i i3 (Ai)7^ i ai+ i(A) which also has a wrong order 
on the rapidities A and Ai. This term as far as / 7i, 2 (Ai)7^ 2ia2+2 (A) |0) are fortunately fixed by 
another set of commutation rules defined by Eq. (14ip as well as by the linear system of equations 
( 14211431) with d\ = and b% — 3. Putting together all the above steps we are finally able to reorder 
the products of creation fields in a suitable manner. 

We now turn to the computations involving Xi+a,a(A)0 2 (A 2 , A 3 ) |0) for d — 1, ■ ■ • , N — a. To 
this end we shall first discuss the cases where d > 3. In this situation the computations are 
somehow simpler since the azimuthal spin of the operator Td+ a ,a(ty exceeds that associated to the 
two-particle state 2 (A 2 ,A 3 ) |0). In fact, by using the commutation rules (I52|I68|I72II126P we can 
reduce the corresponding calculations to the action of several annihilators on the reference state 
|0) and consequently we find, 



The results for d = 1 and d = 2 are obtained after a number of extra steps are performed. 
First we have to commute the operators 7^ + i j0 (A) and 7^ +2a (A) with the fields Ti^(Xj) for 



T d+a , a (A)0 2 (A 2 ,A 3 )|O> = for d > 3, 



(B.4) 
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j = 2,3 present on the first part of the state ^2(^2, A3) |0). This task is accomplished with 
the help of Eqs. (l52|ll26l) which generates the following type of products 7i,i(A2)7^, a+ i(A) |0), 
7i,i(A 2 )7; + i, 0+2 (A) |0), Z J ,a(A) / 7 2i 3(A2) |0) and T a+lia (X)T 2 ,3(X 2 ) |0). We note that these terms pos- 
sess either a diagonal or an annihilation operator that still need to be carried out to the right-hand 
side. This operation for the first two terms is done by using Eq. fl72|) . For the last two terms we have 
to employ Cramer's rule in Eq.(|65l) to compute the commutation rules among the fields 7j, a (A) 
and 7^+i )0 (A) with 72,3 (/i). The last step concerns with the commutations among the operators 
^o+i,a(A), Zj +2 , a (A) and T 13 (X 2 )- This is easily performed with the help of of the commutation 
rules described at end of section [3731 by Eqs. f)68H74")) . Collecting together all the above mentioned 
steps we find that, 

T a+2 ,a(A)0 2 (A2, A 3 ) |0> = Wa+2 (A) Wl (A2) W i(A 3 )o^ a) (A, A 2 , A 3 ) |0> 

+ Wa+1 (A) W 2(A2WA 3 )i^ a) (A, a 2 , A 3 ) ^ 2, ^; : ;^ 3, ^; : ; g(A 2 , a 3 ) io> 

^hA 2 , A 3 J 21 K(A 3 , A2)i]i 
- w a+1 (X)w 1 {X 2 )w 2 {X 3 ) T[ a \X,X 2 ) n ( t +1) {X,X 2 ,X 3 \2) |0) 

+ w a (X)w 2 (X 2 )w 2 (X 3 ) 2 ^ a) (X,X 2 ,X 3 \l)\0), for a<N-2 (B.5) 



T a+1 , a (A)0 2 (A 2 , A 3 ) |0) = T li2 (A 3 ) ^ 1 (a) (A, A 2 ) K +1 (A) Wl (A 2 )P 1 (A 2 , A 3 )P a+1 (A, A 3 ) 

- w a (X)w 2 (X 2 )P 2 (X 2 , A 3 )P a (A, A 3 )] |0) 

+ 7i i2 (A 2 ) [w a+ i(X)w 1 (X 3 )P 1)a+1 (X, A 2 , A 3 ) 

- W a (A)«;2(A3)P2,a(A,A2,A3)]|0) 

+ r a+1 , a+2 (A)w; 1 (A2) W i(A3)o^2 (a) (A, A 2 , A 3 ) |0) 

T a , a+1 (A)o^ 1 (a) (A,A 2 ) [wi(A 2 )™ 2 (A 3 )o?4 a) (A,A 2 ,A 3 |2) 

+ w 2 (X 2 ) Wl (X 3 ) ^ a) (X, ^) ^ 2,h ]li e(X 2 , A 3 ) 



P(A 2 ,A 



3/2,1 



|o> 



+ S 1 a %^ a (X)w 2 (X 2 )w 2 (X 3 ) 2 ^ a 1) (A,A 2 ,A 3 )|0), for a < N - 1 



(B.6) 
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where functions 2 .F 2 (A, A 2 , A 3 |l), Pi, a +i(A, A 2 , A 3 ) and P 2 , a +2(A, A 2 , A 3 ) are given by, 



2< F 2 (a) (A,A 2 ,A 3 |l) 



P(A,A 



\a,3 



\a+2i 0^1 ^(A2> A3) + -P 2 (A 2 , A 3 ) ^ + ^(A, A 2 ) ^ r i ''(A, A 3 



P(A, A 2 )„+2'l 

(a+l)/x ^ \ ^(a) 



o^r +IJ (A,A 2 )o^ 1 w (A,A2)o^i (1) (A 2 ,A 3 ), for a < N - 2 



(B.7) 



P,a+i(A, A 2 , A 3 ) — P a+ i(A, A 3 )o^{ ^(A, A 2 )o^ r { ; (A 2 , A 3 ) — o^ 7 } + j (A, A 3 )o^ r { ''(A, A 2 



r(a+l)/ 



X 



P(A, A 2 )^; 2 P(A, A 2 )°; 2 
P(A, A 2 )a+2|l R(K ^2)1X1,1 



JQ^l ^(A, A3) — iJ 7 ! ; (A 2 , A3) 



(2), 



-R(A,A 2 )^+ 2jl P(A,A 2 )^+ 12 
^(A,A 2 )o+2!i ^(A>A 2 )^i' 2 



a. 3 



P(A, A 2 )o+ 2 !i^(Aj A 2 )o+i|i 



-f2,a(A, A 2 , A3) — P a (A, A 2 )o^ r i ^(A, A 2 )o^ r { j (A 2 , A3) + 0^1 j (A, \3)o^i j (A, A 2 ) 



-(a). 



-(a), 



-^(A, A 2 )°^ 1,1 

-^(A, A 2 )a+i'{ 



^(A,A 2 )f 



^ T0 ^(A,A 3 ), for a<JV-l. 



(B.9) 



R(\ Aa)o+i|i 

It turns out that the above functions can be further simplified as follows. Considering Eqs. 
(lA~T3l EH IA32|) we are able to rewrite 2 ^ a) (A, A 2 , A 3 |l) as, 



-(») M \ \ , , r <.<n,, , x N P(A 2 , A 3 )ij P(A 3 , A 2 ) 2 ;i 



2.1 



2< F 2 (a) (A,A 2 ,A 3 |l) 



3^*2 (A, A 2 , A 3 ) — iJF 2 (A, A 2 ,A 3 



IT#(A 2 , A 3 ) 



P(A 2 , A 3 )^i P(A 3 , A 2 )^ 
+ ^i (a) (A,A 2 )oM a+1) (A,A 2 ,A 3 |2). (B.10) 



By using the Eq. (1A.8IIA.10j) and then the direct comparison between Eq.( 1A.23j) and Eq. (IB. 10ft 
leads us to the identity, 

2 ^ a) (A, A 2 , A 3 |l) = 2 J- 2 (a) (A, A 2 , A 3 ). (B.ll) 

Next, by exploring the permutation property of the two-particle vector 2 (A 2 , A3) (11391) it is 
not difficult to conclude that the right-hand side of Eqs. (1B.5|IB.6j) must be symmetrical under the 
rapidity exchange A 2 <-> A3 and therefore we find that, 



*W(A, A 2 , A 3 ) = 9(\ 2 , A 3 )Pi(A 3 , A 2 )P a+1 (A, A 2 ) ^ 1 (a) (A, A 3 ) 
P 2 , a (A, A 2 , A 3 ) = 0(A 2 , A 3 )P 2 (A 3 , A 2 )P a (A, A^J^A, A 3 ). 



(B.12) 
(B.13) 
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Taking into account the identities flRTTl [RT2llRT3l) we are able to bring Eqs. (1B.5IIB.6l) into 
their more symmetrical form, namely 

T a+2 , a (A)0 2 (A 2 , A 3 ) |0> = Wa+2 (AK(A 2 ) Wl (A 3 )o^ 2 (a) (A, A 2 , A 3 ) |0) 

+ w a+1 (X)w 2 (X 2 )w 1 (X 3 ) 1 ^ a \\, A 2 , A 3 ) ^ 2, ^; : ;^ 3 '^;^ (A 2 , As) |0) 

ft{M, ^3)2,1 K {M, *2)i,i 

(a) /i(As ; A 2 )^'-^ .R(A 2 , A 3 ) 2' 1 

+ w a+ i(A)wi(A 2 )w 2 (A 3 )i^ ; (A, A 3 , A 2 ) ' ' |0) 

tt{*3, ^2>2,1 tt[A2, 

+ ^a(A)u; 2 (A 2 ) W2 (A 3 ) 2 ^ a) (A,A 2 ,A 3 )|0), for a < N - 2 (B.14) 

Ta+l,a(A)0 2 (A 2 ,A 3 )|O) = 



+ 



+ 



+ 
+ 

The commutation of the diagonal field with the last two terms 7i i3 (A;)0i(A.,) and ^^(Ai) have 
been already discussed in the main text of subsection 15.31 Therefore there is no need to repeat the 
procedure used for such two terms once again. Considering the above results together with the 
steps explained for the operators 7i )3 (Aj)0i(A :) ) and 7i i4 (Ai) we find out the action of the operator 



T 1 , 2 (A 3 )o^ 1 (a) (A,A 



, , ,i?(A 3 ,A 2 ) 

1.1 n /\ \ \ 

w a+1 (X)w 1 (X 2 )-^ 7 - x — rTJT P a+1 (\, A 3 ) 



w a (\)w 2 (\2)e(\2, x ^ X " ,h lii p a (\, A 3 



-R(A 3 , A 2 ) 2 ^ 
10) 



r 1 , 2 (A 2 ) ^ 1 (a) (A,A 



-R(A 2 ,A 3 ) 2 ^ 

i?(A 2 , A 3 )-l \ 

w a+1 (\)wi(\ 3 )9(\ 2 , A 3 ) — — r^n p a+i(^ ^2) 



-R(A 2 , A 3 ) 2 ^ 



Wa(A)w 2 (A 3 ) — ^TTjPaiK A 2 ) 



10} 



-R(A 3 , A 2 ) 2 '! 
T a+1 , a+2 (A)^ 1 (A 2 ) Wl (A 3 ) ^ a) (A, A 2 , A 3 ) |0> 



%i,a+i(X)oJ-'[ \x, A 2 )o^ 7 } ''(A, A 



(a), 



-R(A 3 , A 2 )^ 
-R(A 3 , A 2 ) 2 ' x 



wi(A 2 )w 2 (A 3 



i?(A 2 , A 3 )J'J 
w 2 (A 2 )wi(A 3 )— — — 2j^(A 2 , A 3 ) 

^T a _ 1 , a (A)«; 2 (A 2 V 2 (A 3 ) 2 ^ 2 (a - 1) (A,A 2 ,A 3 )|0), for a < N - 1/B.15) 



10) 
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T a ,a(X) on the three particle ansatz |<3> 3 ) is given by 

3 

r a , a {\) |$ 3 ) = w a {\) J] p a (x, \i) |$ 3 ) 



1=1 

-^T a , a+1 (X) J2 E ^(A n , AJ^lA,) f[ ^^^(A^, A^^f '(A, A,) |0) 

l<n<i 2 <3 ii=i fe=l K \ A ih> X h)2,l 

-5X- 1|(l (A) J] £ 2 (A n , A t2 ) W2 (A ;i ) n a.Jo^Ha, A,J |0) 

l<ii<j 2 <3 ( 1= i k=l K y A hi A ik)2,l 

-^- l ' N % :a+2 (X)J2 E 0i(An)^i(A,J^(A J jn ^ 1, ^ ) ^ < (A n ,A ;>fc ) 
x 2 jf(A, A J1; A, 2 ) |0) - ^T a _ lia+1 (A) t 0i(Ah)^i(A j J^(A^|^^ 



R(X\^^ A^)-^ ^ (a— 1) — 1 2 

x6 l <(A il , AjJ— — — -— 2j^<(Ai i; A il )6 l < (A; 1 , Ajji^ (A, A ii; A/J |0) - S a ' Zj_ 2 ,a(A) 
x E E ^i( A ii) w 2(A h )«; 2 (Ai 2 ) JJ ' fc ' 81 ^ fl< (A /fc , A fl ) J r 2 (a ~ 2) (A, A {l , A la ) |0) 

i 1 =lx<j 1 <i 2 <3 k=l K \ X hi X h)2,l 

_5f- 2 ^- 1 ^r a)a+3 (A) Wl (A 1 )«; 1 (A 2 )«; 1 (A 3 )3^ a) (A, X u A 2 , A 3 ) |0) - S^-^T^^X) 

3 2 

x E E u; i( A ii) w i( A i2)^2(Ai 1 ) n <( A 'x»^iJ2^ o ~ 1) (A,Aj 1 ,A i2 ,A ll ) |0) 

i<ji<i2<3 ii=i fc=i 
ii^'i7^i2 



-^> 2 ' JV T a _ 2)a+1 (A) £ £ ^(A J J W2 (A,J^(A^n^<(^'V)^ 2) ( A ' A ^' A ^ A ' 2 ) 1°) 

j 1 = l l<; 1 <( 2 <3 fc=l 

-^ 2 ' 3 T a „ 3 , a (A) W2 (A 1 ) W2 (A 2 )u; 2 (A 3 )o^ 3) (A, A 1; A 2 , A 3 ) |0) , (B.16) 
where functions ^^(X, Ai, A 2 , A 3 ) for c = 0, 1, 2, 3 have been summarized in Eqs fU80H183p . 
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